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DUPLICATE 



Nonlinoar System, Method of Dosign thoroof In the Frequency 
Domain and Computer Program Product 

The present: invention relates to non- linear systems, 
methods of design thereof in the frequency domain and computer 
program products. More particularly, the present invention 
relates to a non- linear system having pre-determinable 
frequency response characteristics. The invention can he 
utilised to design and realise, for example, nonlinear filters 
having a required frequency response or transfer functions 
having specified transfer characteristics or within a control 
system context • 



The possible frequency components in an output signal of 
a linear system are exactly the same as the frequency 
components of a corresponding input signal. Conventional 
linear filter design is based on the principle that energy in 
unwanted frequency bands is attenuated. 



The Dolby filter, which varies the amplitude of the 
output signal as a function of the level and frequency of the 
input, is an example of a nonlinear filter system. However, 
when comparod with the input, the output does not contain any 
additional frequency components. Modulation is another concept 
related to nonlinear filtering, which is associated with 
signal transmission where the signal to be transmitted is 
modulated by a carrier signal and then transmitted through a 
medium. Although a modulation device allows energy to be moved 
from one frequency band to another, the output frequency 
components of such a device depend not only on the input 
components but mainly on the carrier signal. Therefore, the 
energy transfer Implemented by modulation is realised by a two 
input and one output system where one input is the carrier 
signal and the other input is the signal to be processed. 



The prior art lacks a non-linear system and method 
of /apparatus for the design of such a non- linear system for 
predictably transferring energy from one frequency or band of 
frequencies of an input signal to another frequency or band of 
frequencies Independently of any other input signal. Further, 
the prior art lacks a non-linear control system which can 
predictably transfer energy from one frequency band to another 
frequency band. 



It is an object of the present invention to at least 
mitigate some of the problems of the prior art. 



Accordingly, a first aspect of the present invention 
provides a method for designing a non- linear system for 
transferring energy from an input signal having a first 
spectrum at a first pre— determinable frequency or range of 
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freouencies to an output signal having a second spectrum at a 
fSond ^^determinable frequency or range of frequencies. 

Prefera bly ( -the method oomprioee the Gta P s o£ 

identifying the first spectrum of the input signal from which 

energy is to be transferred, 

specifying the second spectrum of the output signal to which 
said energy is to be transferred, and 

calculating, using a frequency domain description of said 
output spectrum expressed in terms of a frequency domain 
description of said input signal and coefficients of a time 
domain description of a generalised non-linear system, the 
coefficients of a said time domain description of said 
generalised non-linear cyetem in order to give effect to the 
energy transfer. 

Advantageously, the present invention 
allows energy at a particular frequency within a given system 
to be transferred to another frequency or band of frequencies 
at which the response of the system is greatly reduced or 
negligible, or 

allows energy of a signal which is transmitted at a particular 
freauency or band of frequencies to be transferred, without 
using an additional modulating signal, to another frequency or 
band of frequencies at which the associated transmission media 
allows signals to pass, or 

allows energy at a particular band of frequencies to he 
transferred and spread over a new wider range of frequencies 
eo as to attenuate the energy by employing the desired 
interkernel and intrakernel effects of nonlinear systems. 



The present invention is based upon the relationship 
between the input and output spectra or frequency components 
of nonlinear systems, and the relationship between the input 
and output frequencies and/or frequency ranges in the 
nonlinear case, in addition, the invention utilises a mapping 
between the time and frequency domains which allows the output 
soectra or frequency content of nonlinear systems to be 
described completely by the coefficients of time domain models 
which represent the filter or non-linear system to be 
constructed . 

A second aspect of the present invention provides a 
method for manufacturing a non-linear system for transferring 
energy from an input signal having a first spectrum at a first 
pre-determinable frequency or range of frequencies to an 

output signal having a second spectrum at a second pre- 

determinable frequency or range of frequencies, said method of 
manufacture comprising the steps of 

(a) designing said non-linear system comprising the steps of 
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identifying the first spectrum of the input signal from which 
energy is to be transferred, 

specifying the second spectrum of the output signal to which 
said energy is to be* transferred, and 

calculating, using a frequency domain description of said 
output spectrum expressed in terms of a frequency domain 
description of said input signal and coefficients of a time 
domain description of a generalised non-linear system, the 
coefficients of a said time domain description of said 
generalised non-linear system in order to give effect to the 
energy transfer, and 

(to) materially producing the non-linear system so designed. 

A third aspect of the present invention provides a data 
processing system which can transfer energy from an input 
signal having a first spectrum at a first pre-determinable 
frequency or range of frequencies to an output signal having a 
second spectrum at a second pre-determinable frequency or 
range of frequencies, said system comprising 

means for identifying the first apeotrum of the input signal 
from which energy is to he transferred, 

means for specifying the second spectrum of the output signal 
to which said energy is to be transferred, and 

means for calculating, using a frequency domain description of 
said output spectrum expressed in terms of a frequency domain 
description of said input signal and coefficients of a time 
domain description of a generalised non- linear system, the 
coefficients of a said time domain description of said 
generalised non-linear system in order to give effect to the 
energy transfer* 



A fourth aspect of the present invention provides a 
computer program product for designing a non- linear system for 
transferring energy from an input signal having a first 
spectrum at a first pre-determinable frequency or range of 
frequencies to an output signal haying a second spectrum at a 
second pre-determinable frequency or range of frequencies, 
said computer program product comprising 

computer program code means for identifying the first spectrum 
of the input signal from which energy is to be transf erred, 

computer program code means for specifying the second spectrum 
of the output signal to which said energy is to be 
trans f erred , and 

computer program code means for calculating, using a frequency 
domain description of said output spectrum expressed in terms 
of a frequency domain description of said input signal and 
coefficients of a time domain description of a generalised 
non-linear system, the coefficients of a said time domain 
description of said generalised non-linear system in order to 
give effect to the energy transfer* 
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«. *^^-k amBC t of the present invention provides a non- 
A fifth a s Pfct ot r ™ energy from an input signal 

i in ?S a y SiSt W snStrSS STSrs^P^determinable frequency 
S V r2nol H "rJKSSS ?o an output signal having a second 
specSSm ll a second pre-determinable frequency or range of 
frequencies, said system comprising 

means for identifying the first spectrum of the input signal 
from which energy is to be transferred, 

means for specifying the second spectrum of the output signal 
to which said energy is to be transferred, and 

for oivino effect to the energy transfer using 
SlSiicSnts of^a time domain description of a generalised 
™« l±«il? svstem said coefficients having been calculated 
^ni S SegSenS? domain description of said output spectrum 
^ooJ S terras of a frequency domain description of said 
ESS 5SnX an™coe1f!c!en?s of'a time domain description of 
a generalised non-linear system. 

Advantageously, the fifth embodiment allows *o processing 
for the determination of the coefficients to be performed of f- 
Sne and merely incorporated into a non-linear system which 
usos th© coefficients. 

Embodiments of the present invention will be described by 
way of examples only, with reference to the accompanying in 
which: 

figure 1 shows the effect of traditional signal 
processing by, for example, a linear filter; 

figure 2 illustrates signal processing according to one 
aspect of the present invention; 

figure 3 depicts a further example of the signal 
processing according to the present invention; 

figure 4 shows a non-linear system arranged to give 
effect to the energy transformation shown in figure 3; 

fiaure 5 illustrates a further energy transformation in 
whioh energy ie distributed over a wider frequency band; 

figure 6 illustrates the power spectral densities for the 
input and output signals of a designed nonlinear system; 

fiaure 7 illustrates the power spectral densities for the 
input and output signals of another designed nonlinear system, 

figure B depicts schematically a non-linear system; 

figure 9 shows a digital implementation of the non- linear 

system shown in figure 8; „«.J4-.i M 

figure 10 illustrates the power spectral densities for 

the input and output signals of the nonlinear system which was 

obtained using a design wnicn improves the filtering effect 

shown in figure 6; 

figure 11 illustrates the power spectral densities for 

the input and output signals of the nonlinear system whioh wee 
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obtained using a design which improves the filtering effect 
shown in figure 7; 

figure 12 ehowe the time domain input and output of the 
nonlinear system with the frequency domain filtering effect 
shown in figure 10; 

figure 13 shows the time domain input and output of the 
.nonlinear system with the frequency domain filtering effect 
shown in figure 11; 

figure 14 depicts the structure of a designed nonlinear 
system; 

figure 15 illustrates the frequency spectrum of a signal 
to be processed using the present invention; 

figure 16 illustrates the result of a Faer Fourier 
Transform of an input signal having the spectrum shown in 
figure 15; 

figure 17 illustrates the results of the n-dimonoional 
(n= 2 and 3) convolution integration for the spectrum shown in 
figure 16; 

figure 18 illustrates the output magnitude frequency 
response of a designed nonlinear system to an input signal 
having the frequency spectrum shown in figure 15; 

figure 19 shows the frequency spectrum of a further 
signal to be processea using the present invention; 

figure 20 illustrates the output magnitude frequency 
response of the same nonlinear system as shown in figure 18 to 
the further signal to be processed having the frequency 
spectrum shown in figure 19; 

figure 21 illustrates the frequency spectrum of a still 
further signal to be processed using the present invention; 

figure 22 illustrates the output magnitude frequency 
response of the same nonlinear system as in figure IB to the 
still further signal to be processed having the frequency 
spectrum shown in figure 21 ; 

figure 23 illustrates the structure of another designed 
nonlinear system; 

figure 24 shows the continuous time realisation of the 
discrete time system in figure 23; 

figure 25 shows a mechanical implementation of the 
continuous time eystem in figure 24; 

figure 26 illustrates the result of a Fast Fourier 
Transform of the further signal shown in figure 19; 

figure 27 illustrates the output magnitude frequency 
response of a further designed nonlinear system to the further 
signal shown in figure 19; 

figure 28 illustrates the output magnitude frequency 
response of a still further designed nonlinear system to the 
signal shown in figure 15; 

figure 29 illustrates the output magnitude frequency 
response of the nonlinear system shown in figure 28 to the 
further signal shown in figure 19; 
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figure 30 illustrates a flow chart for designing a 
nonlinear system according to an embodiment; 

figure 31 illuetrateB a flow chart for designing a 
nonlinear system according to a further embodiment; 

figure 32 illustrates a data processing system upon which 
embodiments of the present invention can be implemented. 

Referring to figure 1, there is shown the principle of, 
for example, traditional low pass, high pass and band pass 
filtering. Figure 1 shows the power spectrum of a signal loo 
both before and after processing. The energy of a signal lOO 
to be filtered comprises two parts, namely, a first part 102 
for further procassing or of interest and a second part 1D4 
which is of no interest. Typically, the second part 104 of 
the signal ic attenuated which raculta in a second signal 106. 
The second signal 106 comprises the original or a copy of the 
firet part 102 and an attenuated portion or an attenuated 
version of the second part 108. 



Figure 2 illustrates the principle of signal processing 
according to one aspect of the present invention. Figure 2 
shows the power spectrum of a signal 200 both before end after 
processing. The signal comprises a first portion 202 and a 
second portion 204. The first portion 202 of the signal 200 
is of interest for further processing or output. Accordingly, 
as a result of the signal processing using the present 
invention, the first portion 202 is retained and the energy in 
the second portion 204 is translated to anotner frequency band 
206. 



T. PET A T I.ED DKCRIPTION 

The theory and method underlying the present invention 
will now be described in general terms in steps (i) to (vi). 

(I) Determine the frequency spectrum of a signal to be 
processed, including the range of frequencies of the signal. 

(II) specify the frequency spectrum of the output signal. 

(ill) Determine the structure of a Nonlinear Auto-Regressive 
model with exogenous inputs (NARX model) to ensure that the 
energy transformation between different frequency bends and 
other design requirements, for example, specifications for 
magnitude and/or phase of the output spectrum over the 
required output frequency band can be met or realised. 

The general expression for a narx model is given by 



6 



□3-03-98 12:30 +01142730312 ' P 15 R-273 Job-835 

O3-03-9B ' 12:29 From:DIBB LUPTfl^LSOP PATENT i TRADE MARKS +011427303J^ T-807 P. 15 Job-082 



0 



y(k> = £y„(iO 0> 
where Y n <*) is a ' NARX nth-order output' given by 

y n (k) = Z S c ..(Wpjll[yik-ii) ffu (k -ij 

P*01i*lp*a =1 1 = 1 i = p + l 

(2) 

with 

p + q = n, l t =l,...,K n , i = l,...,p + q, and s £ " * z] 

li-l 1^=1 

K n is the maximum lag and yM, u(.), and c pq ( J are tho output:, 

input, and model coefficients respectively. A specific 
instance of -the NARX model such 09 

yCk) = 0.3u(lc-l)+0.7y(]c-X)-0.02u(]c-l)u(1t~l)-Q.Q4u(lc-2)u(1c-l) 
-0 . 06y ( Jc-1 )u(lc~3 J -0 • 08y (1c-2 Jy ( k-3 J 

may be obtained from the general form ( 1 ) and ( 2 ) with 

C 0 i<l)=0-3, C 10 (l)-0.7, c 02 <:i,l)=~0.02, c 02 <2,l)— 0*04, 
Cix( l,3)--0.06, C2o(2,3)— 0.08, else Cpg( . )=0 

Simplified designs will be considered below where the 
NARX model with only input nonlinear! ties is employed. 
However, it will be appreciated by one skilled in the art that 
the present invention is not limited to use in relation to 
only input nonlinear! ties. The present invention can equally 
well be used in circumstances of both non-linear outputs and 
non- linear inputs and outputs* Equally, the invention is not 
restricted to realisation as a NARX model. The invention may 
be realised using many alternative model forms either in 
discrete time or continuous time. Models such as the 
Hammerstein and Wiener model, or continuous time modela, for 
example, a nonlinear differential equation model could be used 
or any other model whioh can be mapped into the frequency 
domain. However, for each of the models the main design 
principle is the 'same. 

The NARX model with only input nonlinearities is given by 
equation ( 1 ) where 

^(V-vljflutk-l,) for *>2 

y n (k)=.^ t " Ki (3) 

2 c ioCli>y(k-l 1 )+& )1 (l 1 )u(k-l 1 ) for n=1 
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The structure of the NARX model (l) and (3) is defined by 

the values of W, K n ,n-1 each « (an integer 

between 1 and N inclusive), involves terms of the form 



n 



1=1 



when n > 2 and 

c^UJyCk-lJ , CoidJu^k-lJ, 
when in the model* 



The parameter N in the non- linear model (1) and (3) is 
associated with the realisability of the model required to 
give effect to the energy transformation- The ability to be 
able to realise the energy transformation is determined from 
the relationship between the input and output frequencies or 
frequency ranges of non-linear systems. 



The structure parameters K n , n=l,,~,N, are associated 
with the extent to which specific design requirements such as 
the magnitude and/or phase of the output spectrum over the 
required output frequency hand can he satisfied* Those 
parameters are iteratively determined as part of the design* 
The model could initially be assumed to be of a simple form in 
terras of these parameters. However, if the initial choice of 
parameters does not produce a satisfactory design the 
parameters are progressively or gradually revised according to 
the energy -transfer effect of the resulting non-linear system. 

For systems described by the NARX model (1) and (3), the 
relationship between the input and output frequencies or 
frequency ranges is given by 



(4) 



where fy denotes the range of frequencies of the output, and 
and f Tw _ i denote the ranges of frequencies produced by the 
Nth-order and ( W-l )th-order nonlinear it ies, and 



L)i k 



when 
when 



k«0 



nb 


na 




< 1 


(a + b) 


(a + b) 




nb 


na 




SI 


(a + b) 


(a + b) 






n = N 


and. 



N-l 



where [.] relates to or means take the integer part/ 



8 



03-03-98 12:30 +01142730312 ' P-17 R-273 Job-B35 

03-03-98 ' 12:30 From:DIBB LUPTUMkSOP PATENT t TRADE MARKS +01 U273031JJL T-807 P. 17 Job-082 




i = 



(a+b) 

na-k(a + b) ,nb-k(a + b)] for k = 0,..., i* - 1, 
I.. = [0,nb-i*(a + b) ], 

and the frequencies of the signal to be processed are in the 
range de-fined by the interval [a,b] * 

Given [a„i>] and the required output frequency range £ y , 
the ©max lest N for the NARX model (1) and (3) which can bring 
about the specified frequency domain energy transformation can 
be determined from equation 4, 

(iv) Map the NARX model with the structure given in (ill) into 
the frequency domain to yield the frequency domain 
description. The frequency domain description is given in 
terms of the Generalised Frequency Response Functions (GFRFs), 
H B (jw 1 ,--,jw B ) , n = l,...,N, which, after this mapping, are 
specified in terms of time domain model parameters. 

The mapping of the NARX model ( 1 )and( 3 ) between the time 
domain and the frequency domain is given by 

n_=l,...,N (5) 

Therefore the frequency domain properties of the system can be 
completely defined in terms of the parameters Cpq (.) of the 
time domain description of the system ♦ 

<v) The output frequency response of the non-linoar system (1) 
and (3) is given by 

w 

Y(jw) = £Y n (jw) < 6 > 
where 

Y B (jw) = -^£- jH^jw^.^jwjnutjw,^ (7) 

with 

J(0d<x u 

denoting an integration over the nth-dimensional hyper-plane 
w^,— ,+w n -w 

Based on this relationship, the parameters in 
H n(J w l' - - • * J w n) > n=1 * - ♦ • * W # which, due to the mapping performed 
in (iv), are the same parameters as those in the time domain 
model are determined. This step enables the shape of the 
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output frequency spectrum Y(jw) to be defined which in turn 
ensure a that the epootrum approochea the specified output 
frequency spectrum as closely as possible. 

Different design specifications can lead to different 
implementations of corresponding designs. 

<v*l) Firstly, given knowledge of the input spectrum U(jw) 

and the required output spectrum Y*(Jw). Substituting (5) and 
(7) into (6) yields 

Y(jw) = " ? = T^b^T X, C^d,,— ,!.>* 



> (8) 



n 

f expf-j (Wjl^,— ( +w„l„)]nU(jw i )d<T w 

Equation (8) enables the parameters associated with the 
time domain NARX model, 

c OA fij l n >' ll- 1 K n *«' n=2 w ' 

and 

Ciof^i>r ll-lp * K 1, 

to be determined as follows to implement the required design: 
1 ) Based upon the equation 

. . A l/Vn ^ , . 

n 

J exp[-j (w 1 l 1 + / — ,+w n l a )]n u ^w 1 )dcr w 



w t 4,- v t-w n =u i_1 



determine parameters 

c 0n (il 2 «J ,*n- • - - '^n- 1 ' • • • < K n> • • * ' N * 

using a least squares routine to make the right hand side of 
equation (9) approach the specified output spectrum as closely 
as possible. 

The first term on the right hand side of equation (8) 

1 



jl - i c xo ( l x > exp (-j wl x )J 
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±s omitted from (9). The omitted term represents linear output 
terms in the time domain realisation of the nonlinear system 
and these may not be needed to achieve the design at this 
step, hence they are omitted from (9). 

2) In order to augment the performance of a filter designed as 
above, it may be desirable to also design a suitable linear 
filter H(Jw) to improve the approximation to Y*<j/w) obtained 
in 1) above such that 

n=*l < «) l l ,l B =X 

n 

f ©xp[-j (w 1 i 1 +,.-,+w n i n )]n u (j w i> d < r « 

con achieve a oet*er- approximation -to Y*(Jw)* As pari: of this 
linear design the parameters c 10 (Ii), 1\ *i which were 

omitted from (8) to get (9), can be obtained from the 
parameters of the linear filter. 



Design 1 hereafter illustrates the design of a non-linear 
system using this first case. 

(V.2) Secondly, given the input spectrum U(^w), and a 

specified bound for the magnitude of the required output 
spectrum , Y B * ( w ) . 

A bound Y B (w) for the magnitude of the output spectrum 
Y(jw) of the NARX model (1) and (3) can be expressed as 



according to the result in Billings, S.A- and Lang, Zi-Qiang, 
1996, A bound for the magnitude characteristics of nonlinear 
output frequency response functions. Part 1: Analysis and 
computation. Int. J. Control, Vol* 65, pp309-328, where 

|uh—*|u(jw)| 
* v — * 

n 

denotes the n-dimensional convolution integration for the 
magnitude characteristic of the input spectrum and 

|H n (jw 1# -.-,jW n )£ 

represents a bound for the GFRFs magnitude 
l H B (j w i'— /jw n )| 

with wi, . . • , w n satisfying the constraint wi+, . . . , +w n - W. For 
the NARX model (1) and (3), 
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can be evaluated as follows 



|H B (jw lf -,jw n )|* = 



Y B (w) = 



where 



Combining (10) ana (11) yields 
1 



L-ic^U^expt-jwli) 



i,=i 



•* |u(jw)j 



(12) 



(13) 



Equation (11) enables 

c„ = i:ic 0 »(i^-vi n )i 

and cio Ul>. il =1 *1 *° be determined for shaping the 

bound YB(w) for Y(Jw) in order to make this bound approach 
yB*(w). The procedure which can be used to achieve this is as 
below . 

1 ) Based upon the equation 

use a least squares routine to make the right hand side of the 
above equation approach yB* (w ) especially over frequencies 
or^reqSency ran£e to which the required energy transformation 
needs to be implemented. The coefficients C n n = l,— ,N, in 
equation (13) must be contained to be positive since C „ ±m 
the result of the summation of the modulus of the coefficients 

c 0tl (li l n J. *n 1 n sl K »- 

The first term on the right hand side of equation (12) 
1 



La omitted from (13). The omitted term represents llneax- 
output terms in the time domain realisation. These omitted 
terms may not be needed to achieve the design at this step, 
hence they are omitted from (13). 
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2) If necessary, the approximation to Y B *(w) in 1) above can 
be Qupplemented using a linear filter with a magnitude 

characteristic |H(jw)[ such that 



provides a better approximation to the specified bound and as 
a result Ciq(Ii), Z^-l , • - . , JC^ , which were omitted from (12) to 
get (13), can be obtained from the linear filter parameters. 

Design 2 hereafter illustrates the detailed procedure of 
the above and several examples of the design. 



( V, 3 ) Thirdly, there are many practical situations which 

should be dealt with on an individual basis. There follows 
examples of two such situations » 



(a) Referring to figure 3 there is shown an input signal 
having a spectrum U(Jcd) 300 whlcn comprises a portion 302 of 
interest, between frequencies a and e, and an portion 304, 
between frequencies e and b, to be translated to another 
frequency range 306 defined by frequencies f and 2b while 
retaining the portion of interest 302 within the output signal 
308. 



In order to realise the energy transformation shown in figure 
3, a non-linear system 400 can be constructed as illustrated 
in figure 4 which comprises means 402, 404 and 406 for 
implementing the following components 



Hi(Jw), H 2 Uw l0 jw 2 ), and H(jw). 

h : ( j w ) 402 ana H(jw) 406 can be implemented readily using 
classical linear band pass filters, while Ho(Jw 1 ,j/w2> 404 can 
be constructed in the time domain as shown in the equation 
below 

l^lj-l 1=1 

with the parameters, c 02 (l 1# l 2 , lj = 1,---,K 2 ,1 2 = 1,--,K 2 being 
determined to produce the signal y 2 ( k with a required ' 
frequency characteristic. Each of the components 402, 404 and 
406 have corresponding frequency responses 408, 410 and 412* 
Consequently, the whole non-linear system can be realised as a 
non-linear time domain filter as 

y(k) = H l (q- l )u(k) + H(q- 1 ) £ c^d^l^Ilu (k - l x ) 
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where 



H^q 1 )^ ™ ■ , H (q *)= , 



are the backward shift operator descriptions of the linear 
filters which have the freguency response functions Hi(jw) ana 
H(jw), respectively. 

The equation for y(h) can be further written as 

H lD (q" 1 )H 0 (q- 1 )y(k:) = H 1N (q- 1 )H D (q- 1 )u(k)-hH N (q'' 1 )H lD (q- 1 ) £ C^d^fjudc-l,) 

which is clearly a NARX model that can be described by 
equations ( 1 > and ( 3 ) . 

Although the above shows energy transformation from a 
frequency band to a higher frequency band, energy can equally 
well be transferred from one frequency band to a lower 
frequency band. 

(b) If the objective of the energy transformation and 
hence the desired non-linear system is only to distribute 
energy of the signal to be processed over a wider frequency 
band without energy amplification then a simpler model may be 
sufficient ♦ Referring to figure 5, there is shown 
schematically an input signal 500 In the frequency domain 
comprising energy between frequencies a and b. The desired 
output frequency spectrum 502 comprises two portions, a lower 
frequency portion 504 and a higher frequency portion 506. It 
will be appreciated that the energy of the input signal 500 is 
to be spread over the two frequency ranges 504 and 506. 

A quadratic filter,, for example, 



can be designed to redistribute energy of the original signal 
with frequency components over the frequency band [a,b], 
to the new ranges [0,b-a] and [2a, 2b] without energy 
amplification provided an appropriate a is selected in the 
design prooese described above. 



(vi) The final step in the design process is to 

materially realise, that is, to physically realise the 
designed filter using appropriate software or hardware or 
combination of software and hardware. 
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Although the present: Invention has been described above 
with reference to a narx model, it will be appreciated that 
the present invention is not limited thereto. The method and 
realisation of non-iineer systems having pre-determinabie 
frequency or energy transfer characteristics can equally well 
he utilised using other forms of descriptions of non-linear 
systems or models where there exists a requirement to 
transform or transfer energy at one frequency or band of 
frequencies to energy at another frequency or band of 
frequencies. 



The present invention can be used to realise energy 
transformations over frequencies using a nonlinear system or 
to modify energy transformation of an existing linear or non- 
linear system. 



Furthermore, the present invention can be applied or 
utilised in the field of design and realisation of electronic 
circuits or filters- The energy in the input signal at a 
particular frequency or band of frequencies is transferred to 
doeired frequency hands. Similarly, in mechanical syeteme, 
the addition of non-linear mechanisms could transfer the 
energy of a vibration at an undesirable frequency to some 
other frequency. The present invention may also find 
application in the field of fluid mechanics, for example, in 
the effects of flow around objects (e.g. oil platform legs), 
noise in ducting and pipe flow systems. 



Alternatively, the modifications which are required to be 
effected to a Known linear or non-linear system, for example, 
a mechanical system, in order to bring about a particular 
frequency distribution of energy can be determined using the 
present invention and then the linear or non-linear system can 
be so modified* 



XX, EXAMPLES 



TT.l DESIGN 1 



Design no.i shows an example of how the energy of an 
input signal having predetermined frequency components can be 
transferred to other frequencies. 

Consider digitally filtering an input signal 
u (t) = cost + cos2t, (14) 
wherein the sampling period is T = l/1°° s - 

The first step is to determine the frequency spectrum of 
the signal to be processed. The frequency spectrum of the 
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signal to be processed comprises input: frequencies 
w Al = 1 and w a2 = 2 . 

The second step is to specify the output frequency 
characteristics of the filter. Two different filtering 
problems will be considered for this example* 

The specification for the first problem is to transfer 
the energy in u(t) to the output frequency w ao = 0 and the 
specif icot ion for "the second problem ±© to transfer -the energy 
in u(t) to the output frequency w ao = 4 . 

since tne sampling period is T = l/l00 s, t:ne digital input 
to the filter is 

1 2 

u (k) = cos—— k + cos— — k for k=0,l,„. , (15) 
100 100 

the normalised input frequencies are therefore 

w dl = 1/100 and w dS =3/100, 

and the required normalised output frequencies are 

w <to =0 for the first filtering problem and = 4/100 for the 

second filtering problem. 

For this example / the output frequencies produced by the 
nth-order nonllnearity of a non-linear filter are distributed 
uniformly in 



Nli when - < 1 

j« ' <« + »>) L <* + *>) J (l6) 

Cj I, when 1* 1 j. j. 

> 0 1 (a+b) L (a + b) J 



1 2 

where a= j"oo # ^ ~ To~0 ' ^ x ^ means -to trake -fche integer part: of 
x, and 

=[na- j (a + b) , nb- j (a + b)] for j-0, ...,i*-l , 
I 4 . = [o,nb-i*(a + b)], 

with the difference between any two neighbouring frequencies 
being A = ~^ = ^- 



For n»l, it can be evaluated from the above equations 

that 
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w _r_i 2_l (i7> 

^"Uoo'iooJ 

ana tHe corresponding output; frequencies are 

\ioo'iooj 

For n-2, it can be similarly obtained that 

and the corresponding output frequencies are 

f _1 2 3 4_1 

\°'ioo'ioo'ioo'iooj 

Therefore, a NARX model of up to second order 
nnniinearitv will be sufficient to realise the energy 
??SnsrSrmaSoSs retired by the filtering example thereby 
addressing the third step in the design process. 

According to the above analysis, select a NARX model of 
the form 

y(k)=i;c 10 (k 1 )y(k-k 1 ) + i:i:c 02 (k l ,k 2 )u(k-k x )u(k-k 2 ) (IB) 

as the basic structure for the non-linear filter and take 
K 2 =l for simplicity. The parameters Kj, c 10 (l) ,...,c 10 (K x ) and 
cj a (k lf k,) , ^=0,1 and k 2 =0,l, then need to be determined to 
completely specify the design. 

in order to derive the procedures for determining these 
parameters, consider the frequency domain characteristics of 
the filter model (18). According to J -C.Peyton Jones and 
sVl.Blllings, Recursive Algorithm for Computing the Frequency 
SUe^onse of a Class of Non-linear ° i£ f^^ A ^ u ^°" Mod - la ' 
int! J. Control, 1989, Vol.50, No. 5 1925-1940 the 
generalised frequency reaponae functions of this sixrer 
model are 

H^jwJ^O 



tlCttlki.y exp[-3(w l k 1 + w 2 k 2 )] 

j^i- 2 c "( k i) ex p[-X w t +w *K]J 

H.tjwj,— ,jw.) = 0 for ns3 



H 2 (dw ir jw 2 ) = — (19) 



The output frequency response of the filter (18) under 
the input 

1 2 , ^ A(w t ) iWit , 20) 
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with Kb2, AlwJ^l, w i = KoO' w -i = - w f i = ±l,±2,can bo writtan 

\ Z^tw^Afw^H^w^jw^) for w>0 



Y (jw) = X Y B ( jw) = ^( jw) = j t Wit+W * I=H 

- S A < w i,) A < w i l ) H 2^Wi l /jw li ) for w = 0 

(21) 

where i^i, e{— 2,-l,+l,-t-2} , Y(jw) is related to the output 
spectrum Y ( j w) by 

f2Y(jw) for w>0 

Y<jW)= I Y(jw) for w = 0 



and, similarly, Y B (jw) is related to the nth-order output 
spectrum Y„ ( j w ) by 

( jw) for w >0 

Y„<ji • 



. % /2Yjjw) 
,W) "U«w> 



for w = 0 



This Is a specific form of the general expression for output 
frequency responses of non-linear* systems which are given by 

n 

Y(jw)-X Y n(5w) (22) 

fttl 

where 

Y D (jw) = -^~ r jH n (jw 1 ,... / jw n )nu(jw l )da w (23) 

with 

J(,)dcr w 

donoting an integration over the nth-dimensional hyper-plane 
w x +w n = w and N being the maximum order of the dominant 
system nonlinear! ties. 

Substituting (19) into (21) yields 
?(jw) = ± SA K )A(w ii )pL 



i-f^c^kJexpT-j (w ti + w lt )kj 
l 1 x 1 r i 

— 7, -ZZcta^i^j); 2 A K J A(w i.) ®*H- j +w lj k 2 )l 

1- 2Xo(Ki> expf-^w k,]- 1 ^ z «n- la - 



c 02 (k 1/ k 2 ) f (w,k 1F k 2 ) for w>0 



and 



y ( jw) — h tjw)i:i:c 02 (k 1/ k 2 ) f (W/ ^ ,k g } f or w = o 



(24) 



(25) 
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where 



H(jw) = - (26) 

l-Xc^dc^expf-Jw k,] 



f (w,k 1# k a ) = \ J> (w^) A(w la ) exp[- j (w 4 k x +w ij k 2 )] 



(27) 



2 



y ( jw) = < 



Moreover denoting 

i i 

ZZ c oz< k i' k J> f (w.k^k.) for w>0 

*»=o*,=o (28) 



gives 

Y (jw) = H(jw)Y(jw). wSO (29) 

In view of the fact that: H ( jw) ie the frequency response 
function of a classical linear filter end Y ( jw) is a linear 
function of the filter parameters Coatk^k^), 

k x =0,l and k 2 = 0, 1 , the procedures for determining the 
parameters of the non-linear filter with the given design 
requirements and the structure ( 18 ) are given as : 

(a) From the design requirements, determine the desired 
output frequency characteristic Y # (jw) and choose the 
parameters ^(k,,^), ^-0,1 and k 2 =0,l, appropriately to make 
7(jw) approximate Y*(jw) as well as possible- 

(b) Examine the filtering effect of 
y(lc)=ZZ£o2(Wu(k-k 2 )u(k-k 2 ) # (30) 

lc l= 0k r =D 

where Co 2 (k 1# k 2 ) , ^=0,1 , k 2 =0,l, are the results obtained in 
the procedure (a). If the effect is acceptable, then choose 
H(jw)=l so that the filter parameters 

<= l0 < 1 > = c io< 2 ) =*— ' = Cio< K x ) = 0 ; 

otherwise design the classical linear filter H(jw) to make 
H(jw)Y(jw) satisfy the requirements for the output frequency 
characteristics and obtain the filter parameters 
K 1 , C l0 (l) ,c l0 <2) i^.iO^KJ at the same time. 

Therefore, in order to address the first problem to 
transfer energy from u( t ) , w ml - 1 and w -a =2, to the output 
frequency w ft0 = 0 , take 

Y*(jO) = l, Y*(jw 1 )=0 / w t =^/ 1=1,2,3,4. 
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and for the second problem to transfer energy from u(t), 
w «i = 1 and w a2 = 2 / "to the output frequency w. 0 " 4 , take 

Y*(j^)=l. Y*Cj Wl ) = 0, W^^, 1-1,2,3,0. 



100 



The filter parameters c 02 {k 1# k 2 ) , k,=0,l , k 2 = 0,l, can then be 
determined through the group equations 



[ ? ' (3 I^ , ] = SS c -'" t -^ f ' < I55' ,c " k ' ) 



(31) 



Im 



by the least squares method. In (31) f R (,) and f T ( .) represent 
the real and imaginary parts of f(-)- 



Rewrite (31) as 



Y = X0 
where . 



Y = [Re[Y-(jO)],..vRe[Y^j^)]im[r(jO)],-.. / l m [r(j^)]J 
*=[c„(0, 0) /Co2 (l, 1) ,c 02 <0, l) + c 02 (l, 0)] T 

t'iO, 0, 0) f "(0, 1^1) £'(0,0,1) 

2. ' 2 2 

: : T 

4 



(32) 

(33) 
(34) 



x = 



f * l Ttt'°' 0) > ^do' 1 ' 1 *- f " ( rbV'°' 1) 

f'tO/ 0/0) f'(0, 1,1) f'<0, 0,1) 

2 2 2 

! * * 

■ » ■ 

£,( roV'°' 0) ' f,( r^-' 1 - l »' ^i^' 0 ' 1 '-. 
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2 

2cos 



^ op 



2 c o « 



10 0 
2 

1 0 O 

3 



cos 



10 0 
4 

TOO 



i--a 

2 

3 *1 0 0 

COS 

1 

" 1 00 

cos 



2 

+ c o s 

1 



10 0 

+ t± & M 

2 
10 0 



1 

10 0 



10 0 



0^ 
2 

- 2sin 



a ain 7-7-7- 

T — 

1 2 



2sin 



-sin 



10 0 
2 

10 0 
3 

10 0 

_4 

10 0 



100 
-sin 



10 0 



10 0 

-sin 



10 0 
■sin 
2 

10 0 



10 0 



(35) 



Notice that 0 can be written in the form of (34) due to the 
fact that 

f (w,0,l)=f (W,1,0). (36) 

Therefore, the filter parameters 0 for the first 
filtering problem to transfer energy from u(t), 
w m1 = l and w.2 = 2 , to the output frequency w m0 = 0 can be 
obtained as 

4 =(X T X)- 1 X T Y 1 (37) 
where 



.,0,- -, 0,0,—, o 



(38) 



and the parameters 0 for the second filtering problem to 
transfer energy from u(t), vr ml = 1 and w a2 = 2 , to the output 
frequency w h0 = 4 can be obtained as 



0 2 =(X T xr 1 X 1 Y 2 
where 



Y 2 = 



0,-«-,0 1 0,-/0 



The results of computing (37) and (39) are 



and 

4 = 



c\*(i, i) 



C 2 02 (0, 0) 
C? oa (l, 1) 
_C 2 03 (0, 1)+C J 02(1,0) 
respectively . 



1187 
1187 
-2373.8 

-120 6. 2 
-1206.2 
2413.2 . 



(39) 
(40) 

(41) 
(42) 
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Under -the input of (20), the power spectral densities of 
the input and output of the non-linear filter 
i i 

= £ 1 o2(Q,0)u a (k) + cV(l,l)u 2 (k - 1) + [c 1 O2 (l t 0) + c^sCOa^Ck - l)u(k) 

(43) 

which is initially designed to address the first filtering 
problem above, are shown In figure 6, and the power- spectral 
densities of the input and output of the non-linear filter 

i i 

Y 2 ( *> = H S £2 o 2 (k 1# k 2 > u ( k - k x ) u ( k - k 2 ) 

k 1= 0k 9 =0 

= c 2 02 ( 0 , 0)u 2 (k) + c 2 0 *(l, l)u 2 (k-l ) + [c 2 02 ( 1,0) + c 2 O2 (0 / l)]u{k-l) u (k) 

(44) 

which is initially designed to address the second filtering 
problem above, ore shown in figure 7* 

Further improvement to the performances of the filters 
(43) and (44) may be possible. Therefore, a linear filter, 
H(jw), is designed m order to further improve the filter 
performances . 

To improve the performance of the filter ( 43 ) for the 
first filtering problem, H(jw) is designed to be a fifth 
order low-pass type 1 Chebyshev filter with cut off frequency 
0.5 r ad /sec and 0*5dB of ripple in the pass-band to ensure a 
satisfactory frequency response over the paee band. The result 
is 



Hj(jw) = 



l-a L (2) z' 1 a L (6) z" 



= icr 12 



_ 12 0.0174 + 0.0853z- l +0.1812z" 2 +0.l652z" 3 +0.0933z^ + 0*0159z" 



l + 4.9941z _x - 9.9765z" 2 + 9 . 9648z -3 - 4 . 9766Z" 1 + 0.9942Z" 5 



(45) 

To improve the performance of the filter (44) for the 
second problem, H(jw) Is designed to be a fifth order high- 
pass type 1 Chebyshev filter with cut off frequency 3.9 
rad/sec and 0.5dB of ripple in the pass band for the same 
purpose as in the first problem case . The result is 

w/ . ^b 2 {l) + b 2 (2)z- 1 -h-H-b 2 (6)z' 5 

( jW) — ' 



l-a 2 (2) z" 1 a 2 (6) z 



0.9218 


-4 .6088Z" 1 +9.2175z" 2 


-9.2175Z" 3 +4 . 6088z" 4 


-Q.9218z~* 


1 + 4 


.83812"* -9.3 635z _a H-9 


.0613z° - 4*384 6Z"* +0 


. 848 6z -B 



(46) 

The purpose of the additional linear filter ie to 
attenuate unwanted frequency components in the outputs of the 
filters (43) and (44) for the above two filtering problems 
respectively to make the output of the additional filter 
satisfy the corresponding design requirement. 
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It will bo apprQcioted that the filter parameters 
K x/ c 10 (l) /•-#c l0 (K 1 J , associated with the expression 



|l-Scio( k i) ex P[--^ wk i]| ln H ^w) 



can be obtained as a result of 



dividing the denominator o£ H(jw) by the numerator of H(jw) . 
The specific H(jw) is given by H^jw) and H^Jw) in equations 
(45) and (46) for the two filtering problems respectively. 

The general description for the nonlinear filters 
designed as above is 

" 1 



y(k)^X c 1 o<k 1 )y(k-k 1 ) = 



~ Z c io< q'" 1 jy ( k ) = H ( q' 1 ) y ( k) 



"ZZc^k^k^ufk-k^fk-k^ ZE c oj( Jc 1/ k J )q;^qr , u(k) (47) 

where 



H <q 



=[l"i)c 1 o( k iJq" kl ] . (48) 

and q~ l , q^ 1 and qj 1 denote the backward shift operators. 
Another expression in the time domain for (47) *is 

y(k)=H" I (q- 1 )[ Il^tk^k.lq-V 1 ]u(k) (49) 

Therefore, embodiments of the non-linear filters which have 
boon designed can be realised in a manner as shown in figure 8 
or more specifically as shown figure 9 - 

It will be appreciated that in figure 8 the two 
components 8O0 and 802 are represented by 

k t sOk,aO 

and 

° 2(q } " H (q ] l-a(2)q- *(6>q- 

For -the first: filtering problem in this example 

G l (qr 1 ,«fi 1 )=tic , .{k 1 ,ic l )(i l -V 

with c^Oq,^) 's given by (41) and 

« _ b 1 (l) + b x (2) q" 1 4—4-b 1 (6) q' s 
G * {<1 } 1-^(2) q^-^—a^SJq- 5 
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with b l Ci) 's and a^i) 's given by (45) . 

For* the second, filtering problem in this example 

with c^Oq,)^) 's given by (42) and 

. b*(D + b 2 (2) g-^-^ce) q- s 
^ 2Cq l-a«(2)q^— a a (6)q- 5 

with b 2 (i) f s and a 2 (i) 's given by (4 6) . 



Referring to figures 10 and 11, there is shown the 
filtering effects of the filters designed to address the two 
filtering problems. Figure 10 shows the power spectrum 
densities of the output and input of the nonlinear filter which 
is finally obtained for the first 'filtering problem. Figure 11 
shows the power spectrum densities of the output and input of 
the nonlinear filter which is finally obtained for the second 
filtering problem. The effects of the filtering in the time 
domain of the two filters are shown in figures 12 and 13. All 
the responses indicate that the filters substantially satisfy 
the design requirements. 
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TT ■ ?■ DESIGN 2 

Design No ,2 illustrates the detailed procedure and 
several examples for designing a non-linear system for 
transferring the energy of a signal from a first 
prede terminable frequency or range of froquanoioe -to a second 
predetermlnable frequency or range of frequencies such that 
the output frequency response of the nonlinear system so 
designed is within specified bounds* 

"2.1 Detailed Procedure 

(1) Given u(t), the signal to be processed in the time domain, 
the frequency band [a,<±] to which the energy of u(t) is to 

be transferred, and the user specified bound Y a *(w) for the 
output spectrum Y(jw) over [c,d] for the design* 

(2) Sample the time domain signal u(t) with sampling interval 
T to yield a discrete time series {u(lc)} and perform a Fast 
Fourier Transform (FFT) on the series to compute the 
spectrum U(jw) of u(t) as 



where U d [j (•)] is the result of the FFT operating on (u(k)} 
and M is the length of the data used to perform the FFT. M 
is taken as an even number for convenience* 

(3) Evaluate the range [a,b] of frequencies in u(t) as 
where 1 & is an integer such that 

and l a is an -integer such that 

(4) The relationship between the bound of the output spectrum 
Y B (w) , the coefficients of the NARX model 

y(ic)=iy„(k) 



24 



03-03-98 12:30 +01142730312 1 P-34 

03-03-98 * 12:33 Ffom:DIBB LUPTJfl|L$OP PATENT & TRADE MARKS +011427303 



/ 

R-273 Job-835 

T-807 P. 34/88 Job-082 



y.<*H 



2-2 ^»( i i/—' i .)ri u ^- i ») f ° r 

l a =l 1=1 
Sciod^ytfc-ljJ + ScbitliJutk^lJ for n = l 



and the spectrum U(jw) Is given by 



Y B (w) = 



- 2 cjodi 



1 



y — - — 



|u|*~. *|u (3w)| 



- 2 c io^i 



J oup {-j Wlj) 



where 



K K 



n = N 0 ,— ,N 



are parameters associated with the NARX model parameters 

K n*Conlli'—*lu) * l x = l,---,K n ,—,l n =lr-#K nl for n = N fl ,— ,N , 
|u|*--*|u(jw)| denotes the n-dimensional convolution 

n 

integration for the magnitude |u(jw)|of the spectrum U(jw), 
which is defined by 

|U[*- - Ho ( j w )| = [ " ■ f_~ |o ( j v i )]- |u [ 3 (w-wrr-i-w^) idw^-dw^ 

n 

and n 0 = l when the narx model involves nonlinear terms from 
order 1 to N m 



Based upon this expression, the structure parameters 
N and N 0 and tho NARX model parameters are determined as 
below* 



( i ) Evaluate 




IV-i 








r 








k i=D 



nb 



■to) 



nb 



na 



(a + b) _ 
na 

Ca + b) L (a + b) 



<1 



where £ x ] denotes the integer part of x, 

1 "L (ITbl J +1 

I 4 =[na-i (a + b) ,nb-i (a + b)], for i = 0,---,i* - 1 
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l r = [o,nb-i*(a + b)] 

for n=l,2,.«* until a value of n is reached such that 
part of the specified output frequency range [c,d] falls 
into £f n . This value of n is used as the value of N 0 ♦ 

( 11 ) Evaluate 

for n-2 P 3^«*. until a value of n is reached such that 
the frequency range [c,d] falls completely within the 
corresponding f r ♦ This value of n is taken as the 
value of N . 

(ill) Calculate 

|0|*---*|U(jwl| 

n 

to yield 

|U|*...*|tJ (j2* i/Mm|r for i = 0,...,M/2 

V V ' 

n 

using the algorithm 

[U|*-*|0(3 2ir i/MT)J = Tu[i + (y - l) nl(^) '^i = 0 , f 
(G<0) ,- -,U [n (M-l)J - Convjfu ( 0 ) ,U (M - 1 )].-■• ,[v ( Q ) ,— ,17 <M - 1 

for n = N 0 ,--- / N, where Conv( . ) denotes the convolution 

operation and U(«) and U(.) represent the 
intermediate results of this algorithm. 

(iv) Based on the (i d - i c + 1 ) equatlono 



JcMTl , fdMTl 
d| L i d = round „ and round (x ) 



where i c = roum 

means to take the integer nearest to x, use a least 
squares routine to compute 

C n , n = N 0 ,— ,N, 

under the oonstraint that the result a must be 
positive , and then select the NARX model parameters 

K n , c 0ft (l 1# - v l n ) , 1, = l f ...,K B , — ,1. = 1,— .K. 
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for n = N 0 /---/N 

under the constraints on the summation of the modulus 
o£ the aoeff iaients given by 



(v) If necessary, design a classical linear filter, for 
example, a band pass filter which ideally allows the 
frequency response to be unity over the frequency band 
[c,d] and zero beyond to yield the linear frequency 
characteristic 



1- Xc^di) exp ( 



and therefore determine the parameters associated with 
K t/ c^dj , li-lf'-jK,. Otherwise, all of the parameters 
of c 10 (.) can be taken as zero to yield a model having 
no jreg re aai on terms associated with the output. 

( 5 ) Construct a NARX model as shown in figure ( 14 ) using the 

results obtained in the above (lv) and (v). The nonlinear 
system illustrated in figure 14 comprises a nonlinear part 
1400 and a linear part 1402 # It will be appreciated that 
in figure 14 

XI ,2,2 Three flpecific Enamglea 
Example 1 

This example illustrates a further Implementation or 
design of a nonlinear system following the above detailed 
procedure - 

(1) The signal to be processed is given by 
U tr\-7M 1 sin g t -sin 0 t 

with a = 3 , 3, and M u =1.6. The frequency spectrum 

tf(jw) of the signal is shown in figure 15 indicating that 

the real input frequency range is [1,3.3]. The requirement 
for the design is to transfer energy of the original 
signal to the frequency band [c,d] = [5. 6, 7.6] with the bound 
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on the output spectrum magnitude specified to be Y B *(w) = 1.6 
over this frequency band. 

(2) Sample u(t) with sampling interval T = 0.01 sec to produce 

1 sin ttkT-sin B kT 1 sin 3.3xQ.01k-sin 0.01 k 
u(k) = 2M u - - = 2*1.6x-x 

k = -1999,— , 0,---,2000 

and perform a Fast Fourier Transform (FFT) on this series 
(M-4000) to compute 

u -( d f 0 - ^OiifoO 1 - °<~'> 2000 

and then to yield 

l) = ufj — 1 j = 0 . 01U d [ j-^- 1) 1 = -1999,- vO r v 2000 

v MT / A J 4000x0.01 J d ^ J 4000 J 

the result of which/ in the nonnegative frequency range, 
is shown in figure 16. 

Notice the difference between the real spectrum of 
u(t) in figure 15 and the computed spectrum' in figure 16. 
The differences are due to the errors caused by the FFT 
operation. The design should and will be performed based 
upon the computed spectrum to lead to more practical 
results . 

(3) Evaluate the frequency range [a,b] of u(t) from the 
computed spectrum giving 

a = 0.6283 b = 3.7699 

(4) Design the system structure and parameters 

(i) Determination of N 0 . 

Clearly, the output frequency range contributed by 
the linear part when the input frequencies are within 

[a,b] = [0.6283, 3.7699] is 

f^ ==[a,b]-[0.6283, 3.7699] 

The frequency range £^ produced by the second order 
nonlinearity in this case is obtained as follows. 

Since n — 2, 

nb f na 1 2x3.7699 f 2x0.6283 



f na 1 2x3.7699 T 

\ L (a+b) J~ (3.7699+0.6283) |_ (3.' 



(a+b) L (a+b) J (3.7699+0.6283) |_ (3.7699+0.6283) J 

=1.7143-[ 0.2857 ]=1.7143~0>1 
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and 



L « J_^L_1 +1 J 2x0.6283 "L X -[ 0.2857 ]+l=0+l=l 
|_<a+b)J L 3.7699+0-6283 J 

So 

^ HJi ^U»-i few** (a+b]j jljo^b-ife+bj 

=[na # riJJO^±>-(a+b^=[2x0.fi2a3, 2x3.769^j0, 2x3.7699-(0.6283+3.7699} 

= [1.2566, 7.5398flJ[o, 3.1416] = [0, 7.5398] 

thereby obtained contains part of the specified 
output frequency range [c»d] = [5.6, 7.6]. So, N 0 is 
determined to be N 0 = 2 . 

(ii) Determination of N. 

Evaluating f^-f^U^for n=2 yields 

f*l„=2=f«;lX=[0/ 7.5398]|J[0- 6283 ' 3.7699] = [0, 7.5398] 
To evaluate f,|,, =s = f^ljf^ , calculate ^ first. In this 



case 
nb 



r na "I 3x3.7699 T 3x0.6283 ~| 
1 (a+b) J~ (3.7699+0.6283) |_ (3.7699+0.6283) J 



(a+b) L <a + b) 

-2.5714-[ 0.4286 ]=2.5714-0>1 

i- J_2e_l + 1 J 3*0,6283 ] + !=[ 0.4286 ] + l = 0 + l = l 
L(a*+b)J L 3.7699+0.6283 J 1 

So 

^ -U x i = I oU I i -[na/nbUo^nb-la+bl] 

= [3x0.6283, 3x3.7699jJ[0, 3x3. 7699-{0. 6283 +3.7699)] 
= [1.8848, U.3097jJ[0, 6.9115] = [0, 11.3097] 
Therefore 

^U^,LK=[0* 7.5398|J[0, ll.3097] = [0, 11.3097] 

fy| n= 3 thereby obtained includes the whole specified 

output frequency range [5.6, 7.6]. N is therefore 
determined to be N = 3 . 

(iii) Calculate |U|*—*|u ( jw)| for n = N 0 = 2 and n=N=3, 

> ^ > 

n 

respectively to yield 
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llt^Huf^V^H^Hu^V^O 00 ^-^ i=0/->400qf2r n=2 arri 3 

n n 

The results are shown In figure 17. 

(iv) Based on the (i d -i e +l) equations 

Y a *(7T i/2000x0.01) = 1.6=^C 2 |u|* ju( j K i/2 0 0 0 x 0 . 0 1)| 



(2*r) 



^-7C 3 |u(*|u|*ju(j^ i/20O0x0.0l)| 



with 

i c = r ound[cMT/2jr] ~ round[s .6x4000x0.0 1/2 x *r] = 36 
i d = round[dMT/2^] = round[7 . 6 x4000x 0. 01/2xjr]= 48 

use a least squares routine under the constraint of 
nonnegative solutions to compute C 2 and C 3 , that is, to 
determine C 2 and C 3 , under the constraints of 
C 2 > 0 and C 3 £0, to minimises the following expression 

^1.6-— q 1^^^^2000 x0.01)1-^03 MHU*Mj^i/2000xO.Ol)|J 

The results obtained are 

C 2 ~0 and C 3 =3.8367 

(v) Design, optionally, a linear Butterworth band pass 
filter to attenuate the frequency components beyond 

the frequency band [c, d] = [5.6, 7.6] to yield the linear 
frequency characteristic 



10" 3 (0 • 0986- 0 . 1972q~ a + 0 . 098 6q~ 4 ) 
l-3^Q€33ci- 1 + 5.899 8q' 2 -3. 907 6q" 3 + 0 . 9721q~ 4 l q ~ G 
(5) Construct a MARX model as shown in figure 14 with 

N 0 =2, N==3, K 2 =K 3 ^1, c 02 (l,l)^0, c 03 (l, 1, 1) = 3.8367, 
that is, 

n^j 1,-<1 l„=l i^l lj -11,-1 W 
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=z£ ° 0J (ii^ 2 )ri u < k - i i) + 



ill 



and 



zsz coitii/ifijriu^-ii) 

i=l 

=^(1,1)^^-1)4^(1,1,1^^^ 

r K 10"fo.0986-0.1972q"V0.0986q- 4 ) 
^-Zi^^^^j = l-3.9633q^ l +5.8988q" 2 -3.9O76q- 3 +0.9721q M 

which determines the parameters associated with the NARX model 
parameters K x and c^d^, l^l,*-,!^. 

The output frequency response under the given input is 
shown in figure 18 indicating that the energy has been 

transferred to the specified frequency band [c, d] = [5.6, 7.6] 
with the magnitude of the response below the specified bound 
1.6. 

The frequency response of the above design is examined 
for other input signals below. 



In the first case, consider 

2 ** u [\ (a + b)t: I 

u ( t ) - — r 2cos ~ -cosbt-cosa t 

;r(b-a)t a L 2 J 

where a, b, and M u are defined as above. The frequency 
spectrum of the signal is shown in figure 19. The frequency 
range is clearly the same as that of the signal having the 
spectrum given in figure 15 and the magnitude of the spectrum 
also satisfies the condition that 

|U(jw)|^M u -1.6 

This implies that the frequency response of the designed 
system to this u(t) should in theory also transfer energy into 

the frequency band [c, d] = [5*6, 7.6] with the output magnitude 

frequency response being less than Y B *(w) = 1.6 over this 
frequency band. Figure 20 shows this frequency response and 
indicates that the actual result is consistent with the 
theoretical predictions. 



For the second case, u(t) was taken as a random process 
with the frequency spectrum given in figure 21. The frequency 
spectrum is substantially within the frequency range [l, 3/3] 
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with a magnitude of less than 1.6. Therefore/ the same 
conclusion should apply for- the output magnitude frequency 
response of the designed system to this random input. Figure 
22 shows this response and indicates that the energy is 
transferred to a new frequency band of substantially 

[5.6, 7.6]. Note that the magnitude of the output spectrum over 
this frequency band is well below the specified bound 

Y B *(w) = 1.6. This is because of the effect of the attenuation 
which is due to intrakernel interference of the nonlinear 
mechanism. 

The nonlinear filter which has been designed above can be 
represented by the block diagram shown in figure 23. This 
could be realised electronically by following the approach 
used no realise Design 1 in figure 9. However, for some 
applications the design may need to be realised in continuous 
time, using the bilinear transformation 



where T=0.01 is the sampling interval, q" 1 is the delay 
operator, ana s Is the Laplace Transform operator and 
substituting in the discrete time design in figure 23 provides 
the equivalent continuous time system shown in figure 24* 
Simulating the system in figure 24 produces almost identical 
responses as the discrete time system response of figure 18. 

The system in figure 24 could for example be realised 
mechanically as Illustrated schematically In figure 25 where 
the cubic device is either a material which exhibits a cubic 
response or Is implemented as an actuator which taicos u x (t) as 

input and produces an actuation output u 2 ( t ) which is 
proportional to the cubic power of the input. 

One possible application of this design would be In 
vibration isolation. For example it may be required to 
transfer energy from an input frequency range [1,3*3] to the 
frequency range [5.6,7.6] - The new design could he usod to 
achieve this effect. 

Other much more complicated designs for vibration 
Isolation can be achieved based on the present invention. The 
design procedure would be exactly as described above but the 
realisation would Involve the synthesis of dampers, damping 
materials or actuators with the nonlinear dynamic 
characteristics specified by the designs . 



Btttlirol fl 3l 

Example 2 shows an application of the design above to the 
attenuation of signal energy oyer unwanted frequency bands 
using designed nonlinear effects. 



l + (T/2)s 
l^(T/2)s 
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When using linear structures for the attenuation of 
signal energy in a physical system/ the attenuated energy is 
usually absorbed by devices such as dampers in mechanical 
systems and resistors in electronic circuits and transformed 
to other energy forms such as thermal energy. This may lead to 
undesirable effects and measures , such as using radiating 
devices, sometimes have to be taken to compensate for these 
effects- When a nonlinear system is employed/ instead of 
attenuating the signal energy directly as in the linear case, 
the signal energy at a frequency of interest can be spread 
over a wider frequency band and attenuated by means of the 
counteractions between different terms which compose the 
output spectrum. This means that, to a certain extent, a 
nonlinear design for signal energy attenuation can reduce the 
requirements for using energy absorption devices and is of 
great benefit in practical applications. 

Anotner important application would be, for example, in 
the design of the foundations or the modification of "the 
characteristics of buildings and structures which are in 
earthquake zones* The objective in such an application would 
be to design materials or actuators, tuned as required for 
each structure, which transfer the damaging input energy from 
an earthquake to another more acceptable frequency band or 
spread the energy to be within an acceptable bound over a 
desired frequency range. Spreading the energy using the 
present design should produce significant reductions in 
earthquake damage . 



(1) Design a nonlinear system to attenuate th© energy of the 
signal 

u(t) = ""TTTT 2cos[ -"4r~t j — cosa t-cosb t . 

*r{a-b)t 2 [ v 2 J J 

where a = 3.3, b = l, M u ~1.6. 

(2) The spectrum of the signal Is evaluated. This is 
implemented by sampling the signal with sampling interval 
T = 0.01 sec to produce 

J , 0 -^"< 1 CT>' [ 2C ° a (^ kT )- C ° S '' kT - C ° S * kT ] 

k = -1999,—, 2000 

and perform a Fast Fourier Transform (FFT) on this series, 
The result of the FFT is shown in figure 26. 



(3) Evaluate the range [a,b] of frequencies in the signal based 
on the computed spectrum. The evaluation gives 

a=1.0996 b = 3.1416 
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This is because the computed spectrum indicates 

|u<jw)|>0.05 for w €[l. 0996, 3.1416] 

and 

|u( jw)| <0. 05 for other w 

(4) Assume that, as part of the design, output frequency range 
is [c,d] = [0, 7.3] and the required bound over this frequency 

band is Y fi *(w) = l . 

(5) Following the same steps for the design of the system 
structure and parameters as in example 1 above, 

N 0 , N, and C n/ n = N 0# -- # N were determined as follows. 

(i)N 0 

N 0 i© determined to be N 0 ~ 1 . This is becauae 
= [a, b] = [1.0996, 3.1416]e[0, 7.3]=:[c, d] 
where £y denotes the output frequency band contributed by 

the system linear part, and part of the selected output 
frequency range falls into the linear output frequency 
range . 
<ii> N 

In this case, it can be obtained using 3=1,0990 and 
b=3.1416 that 

f^ =[2.1992, 6.2832jJ[0, 2.0420] 
£*=[0, 9.4248] 

Thus, if the maximum nonlinear order is taken to be 2, 
then the output frequency range of the system is 

£*L2=£*Lk=[0/ 6.2832] 
If the maximum order is taken to be 3, then the output 
frequency range 

Irl.*-^U^=[°* 9.4248] 
Clearly, ^r| n=3 contains the whole selected output frequency 
range [c,d] = [o, 7 . 3] , N is therefore determined to be N^3* 

(iii) C 0 , n = l,2,3. 

C n , n = 1,2/3. are determined to minimize the 
f ol lowing express ion 

t l-C l |u(j2WMT)|"^C 2 M^i^MT)!--^^ MHM*|l<jar ±Jm^ 

under the constraints of C t ^ 0 , i = l,2,3. 

Notice that for this specific example, 
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M - 4000 



fcMTl JO* 4000x0.011 
i c = round — — = round = 0 

[7.6x4000x0*011 
= round = 48 

L 2* J 



TdMT 
i d = round] ■ 



L 2n 

The solution to this minimisation problem is 
C x =0, C 2 =2.1932, C 3 =6.0550 

(6) Select the NARX model parameters 

Kn/Cont 1 !'"*'^)/ li =l#-*vK Bf —,l n -1,— # K ft# for n=N 0 ,— ,N, 

based on the results obtained in (5) as 

K X -K 2 =K 3 =2 
q ) i(-) = 0 / 

q 52 (l # l) = l-1932 r q >2 (l,2) = q 12 (2,l)^0 r ^(2,2) = -!, 

<^ 3 (1, 1,1)=3.0550, ^(2, 2, 2) = -3, and other c^t.) 's are zero 
and consequently construct a NARX model as 

y< k >=2 £ •£ ^( 1 i^*vi n )n u { k - 1 i>-S c o^ 1 i> u < k - 1 i> + 

Ka K t 2 K } K, k, 3 

ll=l l^LL^l 1=1 4 1 =U a =01 3 al 1*1 

=c^(l,l)u*(k-l)+c^(2,2)u 2 ^ 

= l,l932u 2 (k-l)-u 2 <k-2) + 3.0550u 3 (k-l)-3u 3 (k-2) 

Notice that the selected NARX model parameters satisfy 
the relationship 

i,=i 

Z ZMV^H^'M+I^ 2 ^ -|l.l93^+|-a| =2.1932-C, 

K| IS K, 

and the different signs selected for 

c^U, 1) and c 02 (2,2) 

and for 
c 03 (l, 1, l) and c 03 (2,2,2) 



35 



03-03-98 12:30 +01^730312 * P - 45 ^^ R " 273 Job " 835 

12:36 From:DIBB LU^IlsOP PATENT t TRADE MARKS +0114273^^ T-807 P. 45/88 Job-082 



03-03-98 



are such as to give effect to the intra- kernel and inter- 
kernel interferences, which is to attenuate the energy of 
the input signal. 

The frequency domain response of the constructed model to 
the sampled series of the input signal is shown in figure 27. 
It can be seen that excellent energy attenuation has been 
achieved by the designed system* It will be appreciated that 
the input energy in figure 2 6 has been spread over the 
designed frequency band by the nonlinear filter. 

Example 

This example shows another application of the present 
invention to the attenuation of signal energy over unwanted 
trequency bands using designed nonlinear effects. The example 
also illustrates the effect of the same design on a different 
input signal to demonstrate the effectiveness in energy 
attenuation of the designed system in different circumstances. 



(1) Design a nonlinear system to attenuate the energy of the 
signal 

1 sin a t -sin p r 
2n t 



U(t) = 2M u 



with a ^3.3, fi=l, and M u =1-6. 

(2) The spectrum of the signal is evaluated by sampling the 
signal with sampling interval T-0.01 sec to produce 

1 sin okT-ain B kT 1 sin 3. 3 x 0. Oik- Sin 0.01 k 
U(k) ^^ » -2X1.6X-X ^ 

k = -1999,--, 0 f — ,2000 

and then performing a Fast Fourier Transform (FFT) on the 
obtained time series* The result of the FFT is the same 
as that shown in figure 16. 

(3) Evaluating the range [a,b] of frequencies in the signal 
based on the computed spectrum gives 

a = 0.6283 b = 3.7699 

as the computed spectrum indicates 

|u ( j w)| => O. 05 for w€=[0-6283, 3.7 699] 

and 

|U(jw)|<0.05 for other w. 

(4) Assume that/ as part of the design, it is desired that the 
output frequency range is [c,d] = [0, 10 . 3] and the required 
bound over this frequency range is Y B *(w)-l. 
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(5) Following the same steps for the design of the system, 
structure and parameters as in example 1 above, the 
parameters N 0 , N, and C n , n = N 0 ,---,N are determined as 
follows . 

(i) N 0 

N 0 is determined to be N 0 -l due to 
fy t =[a,t>] = [0.6283, 3.7699]e[0, 10.3]-[c, d] , 

which indicates that part of the selected output 
frequency range falls into the linear output frequency 
range . 

(ii) N 

In this example, N is obtained using a-0,6283 and 
b=3 # 7 699 - Theref ore, 
f,, =[0, 7.5398] 
f, 3 ^[0, 11.3097]* 

Hence, if the maximum nonlinear order is taken to be 2, 
the output frequency range of the system is 

fv| n= a=^lX =[°' 7.5398], 

If the maximum order is taken to be 3, the output 
frequency range 

£*U=£*U^=[°' H.3097]. 

As fyl^a contains the whole selected output frequency 
range [c,d] = [0, 7 . 6] , N, for the design, is determined to 
be N=3. 

(iii) C ft , n = l,2,3 

C A , n = l,2,3 are determined in such a manner as to 
minimise 

Jfi-cJu(j2^Mrj-^c, IcI^jariH-^c, |iM*|i<jar 

under the constraints of C^O, 1 = 1/2,3. in this case 
also M=4000 but 

TcmtI To x 4000 x o . oil n 
i c = round[— J = round[ - J = 0 

rdMT"l f lO. 3x4000x0. 01 ~1 

The solution to this minimisation problem is 
C x =0, C 2 -0.2928, C 3 =* 0.9763 

(6) The NARX model parameters 

Kn^CoJl^-vln)/ l a = l,—,K B ,— ,l n =1,— ,K n , for n = N 0 ,--,N, 
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are selected based on the results obtained in (5) as 
K 1 — ■ K 2 = K 3 = 2 

^(1,1) =0.1928, ^(1,2)^-0.1, q^D^O, q H (2,2)=0, 
c^d, 1,1)= 0.6763, c M (l,2,2)=-0.3, and other c^t J »s are zero, 
and consequently a NABX model is constructed 

y<k)= Z S-S ^(ll'-V^)]!"^- 1 !) 

= c 02 (l # l)u 2 (lc-l)+c 02 (l,2)u(k-l)u(k-2) + c 0 ,{l,l,l)u 3 (k-l) 
+Co3(l,2,2)u(k^l)u 2 (k-2). 

=0.1328u 2 (k-l)-0.1u(k-l)u(k-2)+0.6763u 3 (k-l)-0.3u(k-l)u z (k-2) 

The selected NARX model parameters satisfy the 
relationships 

S|c ol (l l )| = 0 = C l 

and the different signs selected for 

c 0? (l,l) and c 02 (l,2) 
and for 

003(1,1,1) and c 03 (l,2,2) 
are in order to give effect to the intra-kernel and inter- 
kernel interferences, which is to attenuate the energy of 
the input signal. 

The frequency response of the constructed model to the 
input signal specified above is shown in figure 2 8 which 
indicates that the required energy attenuation has been 
realised* 

The frequency response of the design above to the input 
signal in example 2 gives the result shown in figure 29. The 
nonlinear system design above clearly also works for the 
signal in example 2 in energy attenuation although the model 
was not specially designed for this signal. This is 
reasonable since the magnitude of the spectrum of the signal 
in example 2 is less than the magnitude of the spectrum of the 
signal in this example over almost all of the input frequency 
band and over the other frequency bands the magnitudes of the 
spectra of the two signals are all zero. This illustrates that 
the above design is effective not only for the input based on 
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which the design is implemented but also for other inputs with 
magnitude frequency characteristics less tfran the magnitude of 
the spectrum of the considered input. 

Ill, FLOWCHARTS 



Referring to figure 30 s there is shown a flow chart 3000 
depicting the steps of an embodiment of the present invention 
in which the output spectrum is specified over- a given range 
of output frequencies. 



The signal to bo prooeeeod and the doeirod frequenoy 
response of the non- linear system to be designed are input via 
etdpfi 3002 to 3006. 

At step 3O02, a digital input signal {u(k)} and its 
sampling period T are to be provided. The range of output 
frequencies [c,d] over which the energy of the input signal is 
to be transformed Is given in step 3004, The output frequency 
range is specified using beginning and end frequencies c and d 
respectively. The distribution of the energy over the output 
frequency range [c,d] is requested or specified in step 3006. 

The frequency characteristics of the input signal are 
determined at steps 3008 and 3010. More particularly, the 
frequency components of the digitised input signal, {u(k) }, are 
calculated using a Fast Fourier Transform at step 3008. The 
range of frequency components contained within the input 
signal is determined from the FFT at step 3010* 



Referring to steps 3012 to 3024, the orders of the 
nonlinearities required to realise a desired nonlinear system 
and, hence, the energy transformation, are calculated. 

Variables n and N 0 are set to one and zero roepoctivdly 
at step 3012 ♦ The output frequency components resulting from 
the nth-order of system nonlinearities are determined at steps 
3014 and 3016. 



A determination is made at step 3018 as to whether or not 
N o -0, that is to say, whether or not the smallest order of 
ay 3 torn nonlinearities which makes a contribution to the 
desired energy transformation has been determined. If N c is 
equal to zero, a determination is made at step 3020 as to 
whether or not part of the specified output frequency range 
falls within fy n * Xf the determination is negative, 
processing continues at step 3026 at which the value of n is 
increased by one. However, if the determination is positive. 
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the value of N Q is set to equal n at step 3022 and processing 
continues at step 3026. 

If N 0 is not equal to zero, then a determination is made 
at step 3024 as to whether or not the specified output 
frequency range lies completely within fy « fy n o fy n -i* If 
the determination is positive, processing continue© with step 
'3027. However, if the determination is negative, the value of 
n is incremented by one at step 3026 and processing continues 
at Step 3014. 

Steps 3028 to 3034 determine the values of the time lags 
of the nonlinear model for all values of n = Nq, Nq + 1'--" N 
and determine the parameters of the nonlinear* system to he 
designed. 

Step 3036 represents a step for fine tuning of the output 
or frequency response of the designed non-linear system. 

Steps 3038 and 3040 determine whether or not the 
resulting frequency response is sufficient to meet the design 
requirements. If the frequency response of the designed 
filter is acceptable, the discrete time version of the filter 
is output at step 3042. However, if the designed non-linear 
system does not have an acceptable frequency response, then 
step 3044 increments the value of K n by one for all n= N(>, 
No+l'*--' N and S * G P S 3030 to 3040 are iteratively repeated. 

The determination in step 3040 as to whether or not the 
designed filter Is acceptable is made, for example, by 
calculating the difference between the required output 
spectrum and the real output spectrum of the designed eyfitem 
over the frequency band [c,d] . If the modulus of the 
difference is below a predeterminable threshold over all the 
output frequency band [c,d], then the designed nonlinear 
system may be deemed to be acceptable. However, if the modulus 
of the difference is greater than the threshold at any 
frequency over [c,d], then the design is refined. 

Referring to figure 31, there is shown a flow chart 3100 
for implementing computer code according to a second 
embodiment of the present invention. 

The signal to be processed and the specification for the 
desired bound to be Imposed upon the output signal frequency 
characteristics is input at steps 3102 to 3106. 



At step 3102, a digitised input signal and its sampling 
period T are to be provided. The range of output frequencies 
[c,d] over which the energy of the input signal is to be 
transformed is specified in step 3104. A bound Y B ( w) for the 
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distribution of the energy over the output frequency band 
[e,d] is input at step 3106. 

The frequency characteristics of the input signal are 
determined at step 3108 and 3110. More particularly, the 
frequency components of the digitised input signal, {uOO}, are 
calculated using a Fast Fourier Transform at step 3108. The 
range of frequency components contained within the input 
signal is determined from the FFT at step 3110. 

Referring to steps 3112 to 3124, the orders of the 
nonlinearities required to realise a desired nonlinear system 
and, hence, the energy transformation, are calculated. 

Variables n and N Q are set to one and zero respectively 
at step 3112. The frequency components resulting from the 
nth-order of system nonlinearities are determined at steps 
3114 and 3116. 

A determination is made at step 3118 as to whether or not 
jj . - O that is to say, whsthar or not the emalieet order of 
system' nonlinearities which contributes to the desired energy 
trans formation has been determined. If N Q is equal to zero, a 
determination is made at step 3120 as to whether or not part 
of the specified output frequency range falls within fy n . If 
the determination is negative, processing continues at step 
3126 in which the value of n ia increased by one. However, If 
the determination is positive, the value of N Q is set to equal 
n at etep 3122 and processing continues at step 3126. 

If, at step 3118, it is determined that N 0 is not equal 
to zero, then a determination is made at step 3124 as to 
whether or not the specified output frequency range lies 
completely within f Y = fy n ^ *Yn-l- if ^ * e ^f mina1:±on i j" 3 
positive, processing continues with step 3127. However, if 
the determination is negative, the value of n is incremented 
by one at step 3126 and processing continues at step 3114. 

Steps 3128 to 3134 determine the values C n for all values 

of n = N 0 . Nn+i N. which represents the summation of the 

modulus of the values of the parameters of the nonlinear 
system . 

At step 3136 the parameters of a NARX model are selected 
given tne constraints imposed upon the non-linear system 
obtained in steps 312B to 3134. A determination is made via 
steps 3138 and 3140 as to whether or not the frequency 
response beyond the frequency band [c,d] of the designed 
filter is acceptable. If the frequency response is acceptable, 
the filter design parameters are output at step 3142. 
However, if the filter characteristics are not acceptable at 
frequencies outside the frequency range [c,d], then a 
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conventional filter, . JKq" 1 ), is designed, at: step 3144, in 
order "to reduce or obviate tho £raquancy responae o£ the 
designed filter outside the range of frequencies tc,d]. 
Finally, the design is completed at etep 3146 by combining the 
designed non-linear and linear filters, if any. 



The determination in etep 3140 aa to whether or not the 
frequency response of the designed filter is acceptable 
outside the range of frequencies [c,d] is made, for example, 
by comparing the modulus of the frequency response beyond the 
frequency range [c,d] with a prede terminable threshold. If the 
modulus is below the threshold over all the frequency range 
outside [c,d], then the designed nonlinear system may be 
deemed to be acceptable. However, if the modulus is greater 
than the threshold at any frequency beyond [c,d], the design 
is refined. 



Referring to figure 32 there is shown a data processing 
system upon which embodiments of the present invention, that 
is to say, the non-linear systems and the methods of design 
thereof can be implemented or realised. It will also he 
appreciated that the non- linear systems so designed may be 
Implemented in digital form using the data processing system 
or other suitable hardware and software. The data processing 
system 3200 comprises a central processing unit 3202 fox- 
processing computer instructions for Implementing the design 
of nonlinear systems given an input signal and particular 
output signal requirements. The data processing system also 
comprises a memory 3204 for storing data to be prooeesed or 
the results of processing as well as computer program 
instructions for processing such data, a system bus 3206, an 
input device 3208, an output device 3210 and a mass storage 
device, for example, a hard disc drive 3212. 



42 



03-03-98 12:30 +01 1427^0312 1 P. 52 R-273 Job-835 

03-03-98* 12:38 F rom:D I BB LUPTONAL SOP PATENT A TRADE MARKS +0114273031^ T-807 P. 52/88 Job-082 



PTONALSOP PATENT & TRADE MARKS +01142730312-. 

# • 



CLAIMS 

1 . A method for designing a non- linear system for* transferring 
energy from an input signal having a first spectrum at a first 
pre-determinable frequency or range of frequencies to an 
output signal having a second spectrum at a second pre- 
determinable frequency or range of frequencies, said method 
comprising the steps of 

identifying the first spectrum of the input signal from which 
energy is to be transfer red , 

specifying the second spectrum of the output signal to which 
said energy is to be transferred, and 

calculating, using a frequency domain description of said 
output spectrum expressed in terms of a frequency domain 
description of said input signal and coefficients of a time 
domain description of a generalised non-linear system, the 
coefficients of the time domain description of the generalised 
non- linear system in order to give effect to the energy 
transfer. 



2« A method as claimed in claim 1 # wherein the frequency 
domain description of the input signal is U(Jw), the time 
domain description of said generalised non-linear system ie 
given by the generalised NARX model 

y(k) = £y n (k) (cl) 

0-1 

where 

y„<k) = £ 2 C pq (l^...,l p , q )]^[ y (k -lj-pf u (lc (C2) 
with 

K„ K„ K n 

p + q = n, 1 4 =l,... f K B , i = l,...,p + q, and £ - 2 

the frequency domain description of the output of the the 
generalised non-linear system is given by 

N 

Y<jw) = 2> B (jw) (c3) 

n-1 

where 

Yn(jW> = 7^^ r J H n<3 w i— ' jwjnu(j Wi )dcr w (C4> 
N is the maximum order of dominant system nonlinearitles, 
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denotes an integration over the nth-dimensional hyper-plane 
w : +,-,+w n = w , and H^w^-" ,jw n ) , n = i,...,N are generalised 
frequency response functions of the the non-linear system. 

3. A method as claimed in claim 2, wherein the mapping from 
the time domain description of the generalised non-linear 
system to the frequency domain description of the system is 
given as 

l-^c l0 (l x ) exp[- j (w^,— ,+w n ) ij W B (jw 1# — ,jw n ) 
ii=i J 

= £co B (l lf — ,1J exp[- jfw.l^.-^+w^J] 

n-1 n-<j Kf>*q \1 / \ 

q=i p=l l l A >mi =l 



(C5) 



P^2 l t A p =l 

where 

H np ( jw^-v jw n ) = ^(jw^-owj H^^t jw 1+1 ,—, jwj expl- j (w^,—,^) l p I 

(C6) 



4. A method as claimed in claim 3, wherfein the relationship 
between the input and output frequencies or frequency ranges 
is given by the following 

where denotes the range of frequencies of the output, and 
fy. and ^ denote the ranges of frequencies produced by the 
Nth-order and ( N -1 )th- order nonlinear i ties, and 



(A 

IK 



when 



when 



nb 



(a + b) 
nb 



na 



L (a + b) J 
na 



<1 



SI 



<C8) 



u«o (a + b) L (a + b) J 

n = N and N-1 

where [.] relates to or means take the integer part/ 

na 



+ 1 



(a + b) J 

I k =[na-k (a + b) ,nb-k(a + b)] for k = 0,„.,i* -1, 
I r = [0,nb-i*(a + b) ], 



and the frequencies of the signal to be processed are in the 
range [a,b]» Given [a,b] and the required output frequency 
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range fy, thejaethod further comprises the step of determining 
the amcixiest N from the relationship above for the generalised 
non-linear system which can bring about the specified 
frequency domain energy transformation. 

5. A method as claimed in claim 4, wherein, with N having been 
determined and K n ,n = 1,...,N , being given a priori, the method 
further comprises the steps of: 

taking N ee Sr and determining the coefficients of the time 
domain model of the generalised non-linear system from the 
description for the system output spectrum given in terms of 
the spectrum of the input signal and the coefficients of the 
time domain model of the said generalised non- linear system* 

6. A method as claimed in claim 5, further comprising the 
steps of 

substituting H n ( jw 1# — , jwj given in (C5) into <C4), and 
substituting the resultant expression for Y n (jw into (C3> to 
obtain the description for the system output spectrum in terms 
of a function of the spectrum of the input signal and the 
coefficients of the time domain model o£ the said generalised 
nonlinear system. 



7- A method for realising or manufacturing a non-linear system 
for transferring energy from an input signal having a first 
spectrum at a first pre -determinable frequency or range of 
frequencies to an output signal having a second spectrum at a 
second pre-determinable frequency or range of frequencies, the 
method comprising the steps of 

(a) designing the non-linear system using the method as 
claimed in any of claims 1 to 6; and 

(b) materially producing the non-linear system so 
designed or using the non-linear system so designed to modify 
materially the transfer* function of an existing linear or non- 
linear system. 



8.. A data processing system for designing a non-linear 
system for transferring energy from an input signal having a 
first spectrum at a first pre -determinable frequency or range 
of frequencies to an output signal having a second spectrum at 
a second pre-determinable frequency or range of frequencies, 
said system comprising 

means for identifying the first spectrum of the input 
signal from which energy is to be transferred , 

means for specifying the second spectrum of the output 
signal to which said energy is to be transferred , and 



means for calculating, using a frequency domain 
description of said output spectrum expressed in terms of a 
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frequency domain description of said input and coefficients of 
a time domain description of a generalised non-Hneor ayotem, 
the coefficients of a said time domain description of said 
generalised non-linear system in order to give effect to the 
energy transfer* 



,9. a system as claimed in claim 8, wherein the frequency 
domain description of the input signal is U(jw), the time 
domain description of the generalised non-linear system is 
given by the generalised NARX model 

N 

where 

y„oo = Z S c M u 1 ,...,i p , q )rty(ic-i 1 ) Ffu(k-i t ) (cio) 

with 

p + q = n, 1 4 =1„..,K B/ i-l,... f p + q, and J) » ]£) 



the frequency aomam description of the output of the the 
generalised non-linear system is given by 

N 

y <jw) = 2> n <jw) (cii) 

where 

Yntjw)^-^— - /HJjw 1# . w jw p )nU(3w 4 )dOL (C12) 

N is the maximum order of dominant system nonlinearitles f 
J( -) dcr w 

denotes an integration over the nth-dimensional hyper-plane 
w l +r- , #+w B =w ^ and H n(3 w i/*"/jw n ) , n = l,— ,N are generalised 
frequency response functions of the said non-linear system. 

10. A system as claimed in claim 9, wherein the mapping from 
the time domain description of the generalised non-linear 
system to the frequency domain desoriptlon of the system is 
given as 
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= £c 0B (l w --vl„)exp[- jfw^^.-^+w.lJ] 

+ZZ S^ 1 !'--' 1 ^^^ 

^1 -~ ■» 1 1 1 » J 



q-1 p-l ij#lp*Q-l 



0 p 

+ Z Z c po( 1 i'—' 1 pK* (:Jw »'" , ^ w » ) 

(C13) 

where 

(C14) 



11. A system as claimed in claim 10, wherein the relationship 
between the Input and output frequencies or fraquoncy ranges 
is given by the following 

f^f^Lk^ (cis) 

where denotes the range of frequencies of the output, and 
and £^ denote the ranges of frequencies produced by the 

Nth-order and ( N -1 )th- order nonlinearities , and 



when 
when 



nb 


na 




<1 


(a + b) 


. (a+b) 




nb 


na 




>1 


(a + b) 


_ (a + b) 






n = N 


and 



(C16) 



where [ . ] relates to or me ana take the integer part/ 



L (a + b) J 



l k =[na-k(a + b) ,nb-k(a + b)] ror k = 0,...,i* - 1, 
I 4 . = [0,nb-i*(a+b) ], 



and the frequencies of the signal to be processed are in the 
range [a,b], given [a,b] and the required output frequency 
range f y, the system further comprises the means for 
determining the smallest N from the relationship above for the 
said generalised non-linear system which can bring about the 
specified frequency domain energy transformation. 



12. A system as claimed in claim 11, wherein, with N having 
been determined and K n? n = 1,.,.,N , being given a priori, the 
system further comprises the means: 
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for taking N as H and for determining the coefficients of the 
time domain model of the generalised non-linear system from 
the description for the system output spectrum given in terms 
of the spectrum of the input signal and the coefficients of 
the time domain model of the generalised non-linear system. 



13. A system as claimed in claim 12, further comprising means 
for substituting H n ( jw l# —, jw n ) given in (C13) into (C12), and 
substituting the resultant expression for Y n (jw into (Cll) to 
obtain the description for the system output spectrum in terras 
of a function of the spectrum of the input signal and the 
coefficients of the time domain model of the generalised 
nonlinear system * 



14. A computer program product for ddcigning a non-linear 
system for transferring energy from an input signal having a 
first spectrujn at a first pre-determinable frequency or range 
of frequencies to an output signal having a second spectrum at 
a aeoond pre-determinable frequency or range of frequencies, 
the said product comprising a computer readable storage medium 
compr i s ing ; 

computer program code means for identifying the first 
spectrum of an input signal from which energy is to be 
trans f err ed , 

computer program code means for specifying the second 
spectrum of an output signal to which said energy is to be 
trans f erred , and 

computer program code means for calculating, using a 
fraquancy domain description of the output epeotrum expressed 
in terms of a frequency domain description of the input and 
coefficients of a time domain description of a generalised 
non-linear system, the coefficients of a the time domain 
description of said generalised non-linear system in order to 
give effect to the energy transfer. 



15. a computer program product as claimed in claim 14, wnerein 
the frequency domain description of the input signal is U(Ji#), 
the time domain description of the generalised non-linear 
system is given by the generalised NARX model 

y(k) = iy„(k) . (C17). 

where 

with 

p + q = n # l^l,...,^, i = l#...,p + q/ and " 5*"" £ 
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-the frequency domain description of the output of the the 
generalised non-linear system is given by 

Y(jw) = 5> n (jw) (C19) 
where 



Yjjw) = 7 J-5r jHJjw, jwjflUfjw.Jda, (C20) 

N is the maximum order of dominant system nonlinearlties, 
J(.) dcr K 

denotes an integration over the nth— dimensional hyper— plane 
w 1 + / --v+w n = w , and H n ( jw 1# -", jw B ) , n = l,...,N are generalised 
frequency response functions of the non-linear system. 

16. A computer program product: as claimed in olaim 15, 
wherein the mapping from the time domain description of the 
generalised non-linoar system -to t:h© frequency domain 
description of the system is given as 



f 



l-Z^ajexpI-j (Wj+,~ ,+w n ) 



x,=i 



, H.(jw l ,...,jw 11 ) 



= Z c ondi/— *l n ) exp[- ^1,+,— ,+w n l n )] 

n-l n-g . . \1 / \ 

+ Z S C po( 1 X^-'lpK^ (jW l'*"' :)W n ) 
p-2 liA p =l 

(C21) 

where 

n-p+l 

H np (jw 1/ »- / jw n )= Z H i^w lf --vjw 1 )H n . l ^ 1 (jw i , 1/ .--,jw n )exp[-j ( Wl 4>— l p J 

( C22 ) 

17* A computer program product as claimed in claim 16, wherein 
the relationship between the input and output frequencies or 
frequency ranges is given by the following 

fr=f*lX. t (C23) 
where fy denotes the range of frequencies of the output, and 
fr_ and fi^ denotes the ranges of frequencies produced by the 

Nth-order and (N-l)th-order nonlinearlties , and 
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i -X 



k=0 



when 
when 



nb 



(a + b) 
* jib 



na 



<a+b) 
na 



<i 



>1 



(C24) 



,ic=o (a + b) L (a + b) . 

n = N and N-l 
where [ . ] relates to or means take the integer part/ 

na 
(a- 

[na-k<a + b) ,nb-k(a + b)] for k = 0,..*,i* -1, 



i — 



na 1 

+1 

+ b) J 



I.. ~[0,nb-i # (a + b) ], 



and the frequencies of tho signal -to bo proooeeed are in -the 
range [a,b], given [a,b] and the required output frequency 
range f Y , the computer program product further comprises _ 
computer program code means for determining the smallest N 
from the relationship above for the said generalised non- 
linear system which can bring about the specified frequency 
domain energy transformation* 

18, A computer program product as claimed in claim 17 , 
wherein, with N having been determined and K n ,n = l,„.,N, being 
given a priori, said product further comprises computer 
program code means of: 

taking N as N and determining the coefficients of the time 
domain model of the generalised non-linear Gystam from the 
description for the system output spectrum given in terms of 
the spectrum of the input signal and the coefficients of the 
time domain model of the said generalised non-linear system. 

19. A computer program product as claimed in claim 18, 
further comprising computer program code means for 
substituting H n ( jw^*--, jw n ) given in (C21) into (C20), and 
substituting the resultant expression for Y n (jw into (C19) to 
obtain the description for the system output spectrum in terms 
a function of the spectrum of the input signal and the 
coefficients of the time domain model of the generalised 
nonlinear system. 



20 ♦ A method for designing a non-linear system for 
transferring energy from an input signal having a first 
spectrum at a first pre-de terminable frequency or range of 
frequencies to an output signal having a second spectrum at a 
second pre-determlnable frequency or range of frequencies 
substantially as described herein with reference to and/or as 
illustrated in the accompanying drawings* 
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21 . a' method for manufacturing comprising the steps of 
designing a non-linear system for transferring energy from an 
input signal having a first spectrum at a first pre- 
determlnable frequency or range of frequencies to an output 
signal having a second spectrum at a second pre-determlnable 
frequency or range of frequencies using a method substantially 
as described herein with reference to and/or as illustrated in 
the accompanying drawings and materially producing the non- 
linear system so designed. 



22 ♦ A data processing system for designing a non-linear system 
for transferring energy from an input signal having a first 
spectrum at a first pre -de terminable frequency or range of 
frequencies to an output signal having a second spectrum at a 
second pre -determinable frequency or range of frequencies 
substantially as described herein with reference to and/or as 
illustrated in the accompanying drawings. 



23. A computer program product comprising computer program 
code means for designing and/or implementing a non-linear 
system for transferring energy from an input signal having a 
first Bpeotrum at a first pre -de terminable frequency or range 
of frequencies to an output signal having a second spectrum at 
a eeoond pre -determinable frequency or range of frequencies 
substantially as described herein with reference to and/or as 
illustrated in the accompanying drawings. 

24. A non-linear system for transferring energy from an input 
signal having a first spectrum at a first predetermlnable 
frequency or range of frequencies to an output signal having 
a second spectrum at a pre-determlnable frequency or range 

frequencies substantially as described herein with 
reference to and/or as illustrated in the accompanying 
drawings . 



25. A non-linear system which can transfer energy from an 
input signal having a first spectrum at a first pre- 
determlnable frequency or range of frequencies to en output 
signal having a second spectrum at a second pre-determinable 
frequency or range of frequencies, said system comprising 

means for identifying the first spectrum of the input 
signal from which energy is to be transferred, 

meane for specifying the second spectrum of the output 
signal to which said energy is to be transferred, and 

means for giving effect to the energy transfer using 
coefficients of a time domain description of a 
generalised non-linear system, said coefficients having 
been oalculated using a frequency domain description of 
said output spectrum expressed in terms of a frequency 
domain description of said input signal and coefficients 
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of a time domain description of a generalised non- linear 
eyetom • 



26. A non-linear system for transferring energy from an input 
signal having a first spectrum at a first pre- determinable 
frequency or range of frequencies to an output signal having 
. a second spectrum at a second pre-* determinable frequency or 
range of frequencies, said system comprising means for 
implementing a method as claimed in any of claims 1 to 7. 

27* A non-linear system for transferring energy from an input 
signal having a first spectrum at a first pre-determinable 
frequency or range of frequencies to an output signal having 
a second spectrum at a second pre-determinable frequency or 
range of frequencies, said non-linear system comprising a 
data processing system as claimed in any of claims 8 to 13* 

28 . An article of manufacture comprising a computer usable 
medium with computer readable program code means embodied in 
the medium for designing a non-linear system for transferring 
energy from an input signal having a first spectrum at a first 
pre-determinable frequency or range of frequencies to an 
output signal having a second spectrum at a second pre- 
determinable frequency or range of frequencies, the computer 
readable program code means in said article comprising: 

computer readable program code means for identifying the 
first spectrum of an input signal from which energy is to be 
transferred, 

computer readable program code means for specifying the 
second spectrum of an output signal to which said energy is to 
be transferred, and 

computer readable program code means for calculating, 
using a frequency domain description of the output spectrum 
expressed in terms of a frequency domain description of the 
input and coefficients of a time domain description of a 
generalised non-linear system, the coefficients of a the time 
domain description of said generalised non-linear system in 
order to give effect to the energy transfer. 
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ABSTRACT 



Nonlinear System, Method of Design thereof in the Frequency 
Domain and Computer Program Product 

The present Invention relates to nonlinear systems and 
methods of design -thereof in t:he frequency domain. Typically , 
conventional linear filter design involves attenuating signals 
at frequonciee whiah are not: of interest and dissipating the 
energy at those frequencies as, for example , heat or sound. 
However, in most systems, it is not always convenient to 
design a linear system or design a system solely with energy 
attenuation in mind. Therefore, tjie present invention 
provides a non-linear system and method of design thereof in 
the frequency domain which can he used to transfer energy at a 
first pre-determinable frequency, or frequency range, to a 
second pare -determinable frequency* or frequency range. Using 
the method of the present invention a non-linear system can he 
developed which can meet given energy transfer requirements or 
a non- linear system can be designed which can alter the 
transfer function of an existing non-linear or linear system. 
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Figure 1 Traditional linear filter designs based on 
attenuating unwanted effects* 



Power spectrum 



2Q0 



^ transfer to another frequency band 




Figure 2 The new design which exploits nonlinear 

properties to transfer the unwanted effects to 
another frequency band 
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Figure 3. A design specification to transfer the unwanted 
frequency components to a different frequency 
band 
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Figure 4 The structure of a non-linear system designed to 
implement the specification shown in figure 3. 
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Figure 5 Illustration of energy distribution over a wider 
frequency band using the present invention 
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Figure 6 The power spectral densities of the output and 
input of the nonlinear filter which is initially 
designed to address the first filtering problem 
in Design 1. 
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Figure 7 The power spectral densities of the output and 
input of the nonlinear filter which is initially 
designed to address the second filtering problem 
in Design 1. 
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Figure 8 The structure of the nonlinear filters designed 
in Design 1 . 
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Figure 9. The block diagram of the nonlinear filters 
designed in Design 1. 



6/26 



' Sp toj 



03-03-98 12:30 +01142730312 ' P 69 R-273 Job-835 

03-03-98 12:41 Frora:OIBB iUPIMLALSOP PATENT 4 TRADE MARKS +011427303^ T-807 P. 69/88 Job-082 



J^^LSOP PATENT 4 TRADE MARKS +01142730^^ 



Pyjy i The Power Spectral Density of Vj 

250 1 ! I i , < i ■ 




i* 1 ' 1 ' i — 1 « - J w 

O 0.6 1 1,5 2 2.5 3 3,5 4 4-5 5 a 

P The Power Spectral Density of u 




Figure 10 The power spectral densities of the output and 
input of the nonlinear riiter eventually 
obtained for the first filtering problem in 
Design l. 
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Figure 11 The power spectral densities of the output and 
input of the nonlinear filter eventually 
obtained for the second filtering problem in 
Design l. 
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Figure 12 Input and output of the nonlinear filter 

designed for the first filtering problem in 
Design 1. 
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Figure 13 Input and output of the nonlinear filter 

designed for the second filtering problem in 
Design 1. 
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Figure 14 The structure of nonlinear systems designed in 
Design 2. 
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Figure 15 The frequency spectrum of the first signal to 
process in example 1, Design 2. 
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Figure 16 The computed result for the spectrum given in 
figure 15, 



\U_ I* - *\U (JO) )1 for n - 2 and 3 




Figure 17 Results of n-dimen$ional convolution 

integration for \U(j(&)\ in example 1, Design 2; 
n=2— solid, n= 3— dashed . 
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Figure 18 Output frequency response of the nonlinear 

system designed In example 1 , Design 2, to the 
first signal with the frequency spectrum shown 
in figure 15. 
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Figure 19. The frequency spectrum of the second signal to 
process in example 1, Design 2. 
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Figure 20 Output frequency response of the nonlinear 

system designed in example 1/ design 2 to the 
second signal with the frequency spectrum shown 
in figure 19. 
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Figure 21 The frequency spectrum of the third signal to 
process in example 1, Design 2. 
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Figure 22 Output frequency response of the nonlinear 

system designed in example 1, Design 2, to the 
third signal with the frequency spectrum shown 
in figure 21. 
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Figure 23 The structure of che nonlinear filter designed 
in Design 2, Example 1 
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Figure 24 The continuous time realization of Design 2, 
Example 1 
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Figure 25 Mechanical realization of Design 2, Example 1 



\U(j<o)\ 




o co (rads/sec) 



Figure 26 The computed result for the frequency spectrum 
of the signal to be processed in example 2, 
Design 2. 
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Figure 27 Output frequency response of the nonlinear 
system designed in example 2, Design 2. 
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Figure 28 Output frequency response of the nonlinear 
system in example 3, Design 2, to the first 
input signal with the same frequency spectrum 
as shown In figure 15. 
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Figure 2 9 Output frequency response of. the nonlinear 

system in example 3, Design 2, to the second 

input signal with the same frequency spectrum 

as shown in figure 19. 
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Figure 30 Plow Chart of* the computor cod© which 

implements the present invention using a 
technique which directly formulates the output 
spectra of nonlinear systems to design. 



3000 



Given a sampled series {u(k)} of 
the signal u(t) to process and 
the sampling period T. 



3002 



(a) 



Given a range lc,cij of frequencies 

specified for the output components of 

the filter to design 



3004 



Given a spectrum Y*(jw) specified for 
the spectrum of the filter output 

over the frequency band [c,d] 



3006 



< 



Calculate U(jw), the spectrum of the input, using 
a Fast Fourier Transformation algorithm to yield 

where Mis the length of tire sampled series ludoT. 



3008 



Evaluate the range [a,b] of frequencies in u(t) from 
the computed input spectrum. 



X 



3010 



Set n=l, N 0 ^ 0 3012 




Perform the following computations 
i =1 1 + 1 



1 [ (a+b) ]" 



J l = [na»i (a + b) ,nb-i (a + b)], 
for i = Q,--*/i*-l 
I.. = [0,nb-i*(a + b)] 
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(c)S 



Evaluate the frequency range fy produced by the 
nth-order nonlinear output as 



i -1 

Ui. 

U=0 



when 
when 



nb 



(a-i-b) 

nb 
(a + b) 



na 



(a+b) 

na 
(a + b) ^ 



< 1 
>1 



where [ x ] denotes the integer part of x. 



n-n+1 



3026 



3016 




Yes 



Whether part of the specified output 
frequency range [c,d] falls into f„ ? .^3020 



NO. 



Yes 



N 0 = n 



3022 



No 



-^Whett 



Whether [c,d] has completely fallen into f y — f y _ Uf Yi 
31 



Yes 
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N=n 
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Set K„=l, for n = N fl> "-, N 3028 



Calculate 

J ey.p[- j (w,l , + w : ,l jjfl u ( j wj do w 

numerically at W - 271 i/MT, i=0,...,M/2, 
for 1, =l,---,K n ,—,l n =!,-••, K n , and n = N 0 ,-.,N 
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(d)< 



(e) 



< 



Compute 

4=r °H>J - L 2n J 

where round(x) means to take th© integer 
nearest to x 



fdMTl 

and i d = round — — I 
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Based on equations 
Y 



( 271 \ A ' l/Vn" A 



i = i c ,i c + 1,— ,i d 

compute 

using a least square routine, , 



3034 



Design a conventional linear 
filter H(q -:l ) to improve the above 
design result within the frequency 
band [c,d] and attenuate the 
responses over the frequency bands 
which are beyond [c,d]. 



Set 


K 


„ =K B +1, 


for 


n 
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Check the frequency domain response 
of the filter 

yCX)-H(q--)£ £ c Dfi (l w »vljflu(Jc-l..> 
over the frequency band [c,d] . 



No 
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Whether the response, in comparisic 
with the specified response Y (jw), 

- — is OK*? 



Yes 



The design is finished and 

the discrete time version of the 

designed nonlinear system is 

y(k) = H(q- l )X £ ^(^^-vljflutk^L) 
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Figure 31 Flow Chart of the computer code which 

implements the present invention using a 
technique associated with the bound for output 
spectra of nonlinear systems. 
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{ Start ) 



Given A sampled acriea {u(k)} of 
the signal u(t] to process 
and the sampling period T . 



(a 



)< 
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Given a range [c,d] of output frequencies 



Given a bound Y (w) specified for the 
spectrum of the output over the 

the frequency band [c,d] 



(b) 



< 
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Calculate U(jw), the spectrum of the input, using 
a Fast Fourier Transformation algorithm to yield 

where Mis the length of the sampled series {u(k)}. 



Evaluate the range [a,b] of frequencies in u(t) from 
the computed input spectrum. 



Set n=l, N n = 0 
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Perform the following computations 
i ~ I — +1 



r na i 



I s = [na-i (a + b) ,nb-i (a + b)l 
for i = 0,---,i'-l 

I r =[o,nb-i*(a + b)] 
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Evaluate the frequency range f». produced by the 
nth-order nonlinear output as 

V-: 

[Jli when 



nb 



when 



1=0 



(a+b) 

nb 
(a + b)' 



na 



L (a + t>) 
na 
(a + b) 



<1 
>1 



where [ x ] denotes the integer part of x* 
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Whether N 0 ^ 0 
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part of the specified output 
frequency range [c,d] falls into f v ? J^12Q 




Yes 



N 0 =» 




No 
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Whether [c,dj has completely fallen into f Y = Uf*^ ? 
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Yes. 



1 



N=n 
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Calculate |u|*- - -* U ( j w)| , the nth-order convolution 



integration on |t)(jw)|, numerically to yield 
|Ul*—*|u(j27i i/MT)|, i = 0,-..,M/2.for n=N Q% ~-iN . 
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Compute 

\ =round — j and i d =round~ 

where round (x) means to take the integer 
nearest to x 



Based on equations 

*-te0-i;5M.M K<^)1 

compute 

C n , n-N or -,N, 
using a least squares routine under the constraint that 
the results must be positive. 



Select NARx model parameters 

K n.C 0n (l l# .-. # l n ) t l x =l,-/K n/ -,l n = l f — / K ||1 

for n = N 0/ ---,N / under the constraints 

and other conatrainta specified by specilic 
design problems. 



Check the frequency domain response 
of the filter 

y(kj=X £ ^(l^ 'Ajfrck-lj 

over the frequency bands beyond [c,d], 
which should, ideally, be negligible. 
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Whether the frequency response 
beyond the band [c,d]is OK? 
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The design is finished and 
the discrete time version of the 
designed nonlinear system is 



yoo=£ i *.(V"A0ftu<k-:L 



n=K 0 1 
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Design a conventional linear band 

pass filter Htq' 1 ) to attenuate the 
responses over the frequency bandfi 

which are beyond [c,d] . 
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Complete the design by constructing 
a nonlinear filter the discrete time 
version of which is 
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Figure 32 The data processing system on which embodiments 
of the present invention can be implemented 



